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Executive  Summary 


During  the  final  perfomance  period,  from  June  14, 1996  to  June  15, 1997 
with  no  cost  extension  to  October  14, 1998,  the  following  accomplishments  were 
obtained.  These  accomplishments  are  delineated  with  respect  to  the  research 
participants  in  the  collaboration. 

The  coupled,  highly  nonlinear  equations  that  govern  the  interaction  of 
plasmas  with  electromagnetic  fields  yields  a  rich  variety  of  physical  phenomenon 
that  are  of  interest  to  applications  in  high  energy,  nuclear  physics,  particle 
accelerators,  high  power  microwave  generators,  and  high  resolution  radar  systems. 
The  theoretical  contributions  that  are  discussed  in  detail  are,  i)  the  stability 
properties  of  coherently  gyrating  relativistic  electron  beams,  ii)  the  chaos  of  halo 
formation  in  the  propagation  of  relativistic  electrons  in  a  periodically  focused, 
solenoidal  channel,  and  iii)  the  thermodynamic  description  of  the  relaxation  of  two 
dimensional  Euler  turbulence  by  using  Tsallis  Statistics. 

The  Vlasov  kinetic  theory  of  the  stability  of  the  radiation  fields  generated  by 
the  interaction  of  a  relativistic  plasma  with  a  solenoidal  equilibrium  field  has 
received  much  study.  These  studies  generally  assume  randomly  isotropic 
distribution  in  the  gyro-phase  angle.  The  work  contributed  by  C.  Chen  and  J. 
Davies  of  MIT,  relaxes  this  assumption  and  considers  gryroscopically  anisotropic 
distribution,  of  temporally  dependent  and  spatially  dependent  form.  The  details  of 
this  work  appears  in  Appendix  I.  The  resulting  coupled  integral  equations,  yield  a 
dispersion  relation  for  the  electrom^etic  modes  that  describes  the  propagation  and 
growth  characteristics  of  the  modes.  Using  the  additionally  simplifying 
assumptions  of  zero  thermal  components  to  the  electron  beam  distribution,  the 
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integral  equations  for  the  jdispersion  xelation  nan  be  xeduced  lo  nlgehraic  equations . 
Specific  equilibrium  distributions  in  the  canonical  coordinates,  yields  a  tenth  order 
polynomial  dispersion  relation  that  is  investigated  numerically.  If  there  is  no  spread 
in  energies  in  the  time-dependent  equilibrium  distribution,  then  the  integral 
equations  reduce  to  three  algebraic  equations  that  relate  the  Fourier  components 
E+(k,w),  E-(k,w-2wc),  Ez(k,  w-wc).  Therefore,  the  eigenmode  system  consists  of  a 
right  handed  polarized  electromagnetic  waves,  a  left  handed  polarized 
electromagnetic  waves  and  electrostatic  waves.  These  components  have  the  same 
wave  number  and  the  same  spatial  and  temporal  growth  rates,  but  they  have 
different  frequencies .  In  the  jcase  of  no  axial  spread  in  momentum,  the  spatial 
dependent  equilibrium  yields  the  following  results  for  the  dispersion  relation.  The 
three  coupled  eigenmodes  are,  E-(k,w),  E4-(k+mw/j)z,w),  Ez(k+mw/pz,w).  Note 
that  the  modes  occur  at  the  same  frequency,  yet  have  different  wave  numbers.  This 
research  continues,  with  the  aim  of  analyzing  the  general  integral  dispersion 
relation. 

The  propagation  of  relativistic  electron  beams  in  a  periodic  solenoidal 
channel  can  result  in  chaotic  particle  orbits,  even  if  the  beam  is  rms  matched  to  the 
propagation  medium.  This  occurs  when  the  space  charge  waves  are  nonlinearly 
coupled  to  the  self  consistently  generated  space  charge  field  of  the  electron  beam. 
For  beams  that  are  space  charge  dominated  in  their  propagation  characteristics,  this 
nonlinear  coupling  results  in  unstable  growth  of  particle  orbits  that  may  be 
considered  chaotic.  These  unstable  particles  occupy  a  halo  region  surrounding  the 
beam  and  represent  increased  emittance  and  radial  transport  as  compared  to  the 
non-space  charge  dominated  beam.  This  work  was  investigated  by  C.  Chen,  W. 
Marable  and  Y.  Fink.  The  result  of  this  numerical  investigation  by  lest  particle 
simulations,  concludes  that  beam  distributions  that  are  peaked  on  axis,  result  in 
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jchaotic  orbits  that ^  jattribntahle  lo m  nnstahlfi  lixed .point  on  Ihe  .b£am  .axis. 
Similarly,  a  beam  distribution  that  is  peaked  off  axis  can  yield  chaotic  orbits  that  are 
attributable  to  unstable  fixed  points  on  the  beam  axis,  and  yields  a  stochastic  layer 
that  outside  the  original  beam  radius.  This  halo  formation  can  be  as  much  as  1.8 
times  the  rms  radius  of  the  beam.  The  details  of  this  investigation  are  summarized 
in  the  publication,  entitled  “  Halo  formation  and  chaos  in  root  mean  squared 
matched  beams  propagating  through  a  periodic  solenoidal  focusing  channel,  by  Y. 
Fink,  C.  Chen  and  W.  Marable.  Found  in  volume  55,  Number  6  of  Physical  Review 
E,  June  1997,  pages  7557-7564. 

The  orbits  from  .a  test  particle  simulation  can  be  .analyzed  to  discern  whether 
the  underlying  dynamics  can  be  characterized  as  chaotic.  Work  performed  by  W. 
Marable  and  M.  Lee,  investigated  the  quantifying  properties  of  dynamical  systems. 
Mr.  Lee,  a  graduate  student  at  Hampton,  successfully  defended  his  Masters’  thesis 
entitled  “  Quantifying  Properties  of  Dynamical  Systems”.  A  product  of  this 
research  was  the  development  and  refinement  of  computer  algorithms  to  extract 
fractal  dimension  and  Lyapunov  spectra  from  time  series  data  sets.  A  summary  of 
these  programs  is  found  in  Appendix  H. 

The  final  research  topic  to  be  summarized  is  the  study  of  Iwo-riimensional 
Euler  turbulence  with  Tsallis  Statistics.  This  work  was  conducted  by  J.  Wurtele 
and  B.  Boghosian.  The  applications  of  these  results  are  of  import  to  both  the  fields 
of  ideal  fluid  flow,  and  to  pure  electron  plasma  studies.  This  work  examines  and 
rectifies  the  experimental  observations  that  two-dimensional  Euler  turbulence  can 
relax  to  a  state  that  does  not  maximize  Boltzmann  entropy,  but  rather  seems  to 
minimize  enstrophy.  In  the  context  of  the  Tsallis  axioms,  a  statistical  theory  is 
developed  and  investigated  numerically.  The  results  show  compatibility  with  the 
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E£&tricted  Minimum  Ensttohpy  jiieory  of  Huang  .and  Driscoll,  nnd  .broadens  Ihe 
application  of  Tsallis  statistics  to  new  fields. 
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Stability  Properties  of  Coherently  Gyrating  Relativistic  Electron  Beams 

J.  A.  Davies"^  and  C.  Chen 
Plasma  Science  and  Fusion  Center 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

ABSTRACT 

An  investigation  is  made  of  the  effects  of  nongyrotropic  equilibrium  distributions  in 
the  phase  angle  of  px  on  the  stability  properties  of  a  relativistic  electron  beam  propagating 
along  and  gyrating  about  an  applied  uniform  magnetic  field.  Perturbations  are  assumed 
to  vary  spatially  only  in  the  direction  of  the  applied  magnetic  field,  so  that  generated 
electromagnetic  and  longitudinal  electric  fields  propagate  parallel  to  the  applied  field.  The 
two  equilibrium  distributions  considered  are  the  time-dependent  distribution  /o{pj.,Pi, 
with  i  =  4>  —  Qct/7  and  the  axial-dependent  distribution  /o(p±,Pz,C)  with  ^  ^  — 

mClcz/pz.  A  Vlasov-Maxwell  analysis  leads  to  integral  equations  relating  the  field  Fourier 
components.  These  equations  reduce  to  algebraic  equations  when  no  spread  in  7  is  present 
in  the  time-dependent  equiUbrium  distribution  and  when  no  spread  in  p^  is  present  in  the 
axial-dependent  distribution.  Numerical  computations  for  these  special  cases  show  that  a 
rich  variety  of  stability  properties  are  obtained  by  changing  the  distributions  in  ^  and 


^  Permanent  Address:  Department  of  Physics,  Clark  University,  Worcester,  Mas¬ 
sachusetts  01610 


AI-1 


I.  Introduction 


During  the  psst  two  decades,  extensive  studies  have  been  made  of  the  stability  prop¬ 
erties  of  a  relativistic  electron  beam  propagating  along  an  applied  uniform  magnetic  field 
If  the  beam  possesses  a  population  inversion  associated  with  the  component 
of  momentum  perpendicular  to  the  field  (pj.),  then  this  system  may  be  subject  to  the 
whistler  and  cyclotron-resonance  maser  instabilities,  which  initiate  the  process  of  convert¬ 
ing  electron-beam  kinetic  energy  into  coherent  electromagnetic  radiation.  The  cyclotron- 
resonance  maser  instability  provides  the  basis  for  existing  and  proposed  electronic  am¬ 
plifiers  and  oscillators  such  as  the  gyrotron  and  the  cyclotron  autoresonance  maser.^^ 
With  the  inclusion  of  thermal  background  electrons  in  addition  to  the  beam,  this  insta¬ 
bility  is  of  interest  in  such  problems  in  space-  and  astrophysics^®”^®  as  the  generation  of 
auroral  kilometric  and  Jovian  decametric  radiation. 

In  this  paper,  we  carry  out  a  stability  analysis  of  this  system  (exclusive  of  an  ambient 
thermal  background)  using  Vlasov-Maxwell  theory.  We  make  the  assumption  that  the 
spatial  variation  of  all  quantities  is  in  the  ^-direction  only;  however,  the  electron  beam  is  not 
assumed  to  be  gyrotropic  in  the  phase  angle  (j)  of  the  component  of  the  particle  momentum 
normal  to  the  2;-axis.  Most  previous  analyses  of  this  system  assume  that  the  equilibrium 
distribution  is  of  the  form  fo{p±,Pz),  i-c.  that  the  beam  is  gyrotropic.  Analyses  of  the 
nongyrotropic  case  are  limited  in  number.  Using  the  eikonal  approximation,  Fruchtman 
and  Friedland^’®  have  considered  the  case  of  a  stationary  amplifier  with  a  nongyrotropic 
equilibrium  distribution  of  the  form  /o(pj.)Pi,C)  —  P'^^ipi.  ~  P±q)^{Pz  ~  Pzo)g{C)^  where 
C  =  <f>  —  mriczfpz  and  flc  is  the  nonrelativistic  cyclotron  frequency.  Kho,  et  aL,^  using 
the  same  equilibrium  distribution,  have  dropped  the  assumption  of  a  stationary  amplifier 
but  assume  that  the  left-hand  polarized  perturbed  radiation  field  can  be  ignored.  Chen,  et 
o/.^®  consider  TE  modes  for  the  case  of  a  helical  relativistic  electron  beam  in  a  cylindrical 
waveguide.  In  an  application  to  space  plasmas,  FVeund,  et  consider  the  case  of  a 
diffuse  electron  beam  propagating  in  a  cold  magnetized  ambient  plasma.  The  beam  is 
coherent  in  phase,  and  has  a  thermal  spread  in  pj_  but  no  spread  in  p^. 

Two  equilibrium  distributions  are  introduced  in  Sec.  II.  These  are  the  time-dependent 
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distribution  fo{p±,Pz.O  where  ^  =  <l>  -  ^ct/lf  and  the  spatial-dependent  distribution 
fo(p±,Pz,  0  where  C  =  By  manipulating  Fourier  transforms  of  the  Vlasov- 

Maxwell  equations,  we  obtain  sets  of  equations  relating  components  of  the  perturbed 
right-  and  left-hand  polarized  electromagnetic  and  the  electrostatic  fields.  For  the  tune- 
dependent  equilibrium,  these  relations  are  given  in  Eqs.  (26)-(28).  The  relations  for  the 
spatial-dependent  equilibrium  distribution  appear  in  Eqs.  (41)-(43).  In  either  case,  these 
relations  are  integral  equations,  not  algebraic  equations.  Eqs.  (26)- (28)  reduce  algebraic 
equations  only  if  there  is  no  spread  in  p  in  7)  iii  the  time-dependent  equilibrium 

distribution.  Moreover,  Eqs.  (41)-(43)  reduce  to  algebraic  equations  only  if  there  is  no 
spread  in  Pz  in  the  spatial-dependent  equilibrium  distribution. 

The  analysis  of  these  integral  equations  is  the  subject  of  current  research  and  is  not 
dealt  with  further  in  this  paper.  Instead,  in  the  remainder  of  the  paper,  we  deal  with  cases 
in  which  the  integral  equations  reduce  to  algebraic  equations. 

Stability  properties  for  the  case  of  the  time-dependent  equilibrium  distribution  with 
no  spread  in  p  are  considered  in  Sec.  III.  Most  generally,  spreads  in  the  pitch  angle 
a  =  (px/Pz)  and  the  phase  angle  <j)  may  be  present.  In  this  case,  the  integral  equa¬ 

tions  (26)-(28)  reduce  to  just  three  algebraic  relations  presented  in  Eq.  (58)  of  Sec.  III.A. 
The  corresponding  exact  dispersion  relation,  relating  complex  frequencies  and  complex 
wave  numbers,  is  the  three  by  three  determinant  relation  in  Eq.  (64).  In  Sec.  III.B,  we 
consider  the  more  restrictive  case  of  no  spread  in  the  pitch  angle  a  in  the  equilibrium 
distribution.  (Then,  both  px  and  have  definite  equilibrium  values.)  In  this  case,  the 
dispersion  relation  in  Eq.  (64)  reduces  to  the  tenth-degree  polynomial  relation  in  Eq.  (69). 
Numerical  computations  of  growth-rate  curves  (Imu>  vs.  k,  real)  and  properties  of  eigen- 
modes  of  Eq.  (58)  are  presented  in  Sec.  III.C  for  the  case  of  definite  equilibrium  px  and 
p^  and  various  equilibrium  distributions  in  0. 

An  analogous  treatment  is  given  in  Sec.  IV  for  the  case  of  the  axial-dependent  equi¬ 
librium  distribution  with  no  spread  in  pz.  Spreads  in  px  and  0  are  still  permitted.  For 
this  case,  it  is  shown  in  Sec.  IV. A  that  the  integral  equations  (41)- (43)  reduce  to  the  three 
algebraic  relations  in  Eq.  (89).  The  exact  dispersion  relation  for  the  system  is  given  by 
the  three  by  three  determinant  equation  (95).  In  Sec.  IV.B,  the  additional  condition  that 
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there  be  no  equaibrium  spread  in  p±  is  imposed.  Then,  (as  in  the  time-dependent  case)  the 
dispersion  relation  reduces  to  the  tenth-degree  polynomial  equation  (100)  relating  the  com¬ 
plex  frequency  and  complex- wave  number.  Numerical  computations  of  growth-rate  curves 
and  properties  of  the  eigenmodes  of  Eq.  (89)  for  this  case  are  presented  in  Sec.  IV.C. 

A  summary  of  our  results  and  conclusions  is  presented  in  Sec.  V. 
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II.  Formulation  of  the  Problem 


A  beam  consists  of  relativistic  electrons,  which  propagate  along  and  gyrate  about  a 
uniform  magnetic  field  Bq  =  Bq&x.  Initially  the  beam  is  in  an  equilibrium  state  in  which 
temporally  and  spatially  varying  electromagnetic  fields  are  absent.  The  initial  growth  rates 
of  these  fields  are  obtained  by  regarding  them  as  small  perturbations  on  the  equilibrium. 
The  system  is  treated  zis  one-dimensional  in  the  sense  that  the  spatial  variation  of  all 
variables  is  in  the  ^-direction  only.  Consequently,  only  electromagnetic  waves  propagating 
parallel  or  anti-parallel  to  the  uniform  field  Bq  are  included  in  this  analysis.  Furthermore, 
equilibrium  self  fields  are  assumed  to  be  neghgibly  small,  so  that  results  discussed  below 
are  limited  to  the  case  of  a  small  ratio  of  the  plasma  frequency  to  the  cyclotron  frequency.^ 

A.  Equilibrium  Distributions 

A  single- particle  momentum  p  can  be  described  by  the  components  pj.)  Pzj  and  (j>, 
where  <j)  is  the  phase  angle  (as  shown  in  Fig.  1).  This  paper  deals  with  two  systems  in 
which  the  phase  angle  <j)  is  not  necessarily  random  in  the  equilibrium  distribution.  One 
simple  constant  of  the  single-particle  motion  involving  ^  is  ^  =  (j)  —  where  Clc  = 
eBo/mc  is  the  nonrelativistic  electron  cyclotron  frequency,  — e  and  m  are  the  electron 
charge  and  rest  mass,  respectively,  c  is  the  speed  of  light  in  vacuo,  t  is  the  time,  and 
1  —  {l  =  (l  +p\/m^c^  +  is  the  relativistic  mass  factor  of  the 

electron.  An  equilibrium  distribution  of  the  form 

/o(p,t)  =  /o(pj.,P^,0.  (1) 

where 

^  =  (2) 

corresponds  to  an  equilibrium  electron  beam  that  is  homogenous  in  the  configuration  space 
a.t  any  given  time.  Since  the  equilibrium  distribution  (1)  is  nonstationary  in  the  momentum 
space,  we  refer  to  it  eis  the  time- dependent  equilibrium  distribution. 
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Another  simple  constant  of  the  single-particle  motion  involving  is 

^  =  - =  (^-m — 2,  (3) 

7  Pz  ^  ^ 

where  v  is  the  electron  velocity.  Use  of  C  gives  rise  to  an  alternative  equilibrium  distribution 


/o(^,p)==/o(Px,Pz,C)-  (4) 

In  this  case,  the  equilibrium  distribution  is  constant  in  time  at  any  given  z.  It  is  analo¬ 
gous  to  a  typical  laboratory  situation  where  the  beam  is  introduced  into  the  interaction 
region  at  some  initial  z  =  0  with  a  given  distribution  in  4>.  Then,  if  interactions  with  the 
electromagnetic  field  are  excluded,  the  electrons  move  along  the  field  lines  with  constant 
Vz  and  gyrate  around  them  with  the  constant  relativistic  cyclotron  frequency  Clf./'y.  Con¬ 
sequently,  the  phase  at  any  value  of  2  is  the  phase  at  2  =  0  plus  We  refer  to  this 

distribution  as  the  axial-dependent  equilibrium  distribution. 

The  time-  and  axial-  dependent  distributions  are  illustrated  with  examples  in  Fig.  2. 
It  is  shown  in  the  Appendix  that  the  two  distributions  in  Eqs.  (1)  and  (4)  are  physi¬ 
cally  different  in  the  sense  that  neither  can  be  transformed  into  the  other  by  a  Lorentz 
transformation. 

In  the  remainder  of  this  section,  we  derive  linearized  equations  relating  components 
of  the  perturbed  electromagnetic  fields  for  each  of  the  distributions  in  Eqs.  (1)  and  (4). 

B.  Perturbation  Analysis  for  the  Time-Dependent  Equilibrium  Distribution. 

The  equilibrium  distribution  is  of  the  form  in  Eq.  (1),  ie.,  fo  (p,t)  =  /q  {p±^Pz,^)  = 
/o  {p±,pz,  —  f^ct/7)  •  The  distribution  is  assumed  to  be  periodic  in  ^  and  (for  fixed  t)  in  <f> 

with  period  27r.  Consequently,  the  normalization  of  /o  {p±,Pz,  <l>  -  0,ct/i)  over  momentum 
space  is  time-independent  and  is  defined  by 

d#±/o  {p±,Pz<,  <l>  -  =  1.  (5) 

As  time  progresses,  the  distribution  evolves  under  perturbations  into 


/  (2,  P,  t)  =  fo  (p,  t)  -f  /i  {z,  p,  t) , 


(6) 
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where  fi  {z,  p,  t)  is  considered  to  be  a  small  perturbation.  To  first  order  in  the  perturbation, 
the  time  evolution  of  this  quantity  is  governed  by  the  linearized  Vlasov  equation, 

=ifi  +  (z,p,t)  -  X  Bo  •  Vp/i  iz,p,t) 

dt  \o  ot  oz  c 

=  e  ^Ei  +  ^  ^C~)  iPx,Py,Pz,t)  ,  (7) 

where  the  total  time  derivative  is  along  a  cheiracteristic  of  the  unperturbed  motion  (ie., 
z  (t  +  r)  =  z  +  v,T,  p,{t  +  T)=p,,px{t  +  r)=Px,4>{t  +  r)  =  <P  +  Qcr /gamma,  and  the 
variable  of  integration  r  is  in  the  range  from  — oo  to  0).  The  fields  Ej.  iz,t)  and  Bi  (z,  t) 
are  regarded  as  small  perturbations  governed  by  Maxwell’s  equations,  which  in  the  present 
treatment  reduce  to 


(8) 

(9) 

^^(z,t)  =47r/)i  {z,t), 
oz 

(10) 

Biz  {z,t)  =  0. 

(11) 

In  the  above  equations; 


Bi±  {z,  t)  —  Eix  {z,  t)  i  iEiy  {z,  t) , 

(12) 

Bi±  {z,  t)  =  Bix  {z,  t)  ±  iBiy  {z,  t) , 

(13) 

Ji±  {z,t)  =  -eno  j  d^pfi  {z,V,t) 

(14) 

Pi  (z,t)  =  -eno  J  d^Pfi  , 

(15) 

p±  =  7mu±  =  p±  exp  {±i<j>) , 

(16) 

where  no  is  the  mean  electron  number  density.  The  fields  (Ei_,J3i_)  and  {Ei+,Bi+) 
represent  right-  and  left-hand  circularly  polarized  transverse  electromagnetic  waves,  re¬ 
spectively,  whereas  the  field  {Eizi  Biz  —  0)  describes  longitudinal  (electrostatic)  waves. 

Because  most  of  the  integration  of  Eq.  (7)  along  characteristics  is  standard,  not  all  of 
the  details  will  be  given  here.  We  remark  that  before  integrating,  it  is  convenient  to  express 
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the  X-  and  y-components  of  the  vectors  that  appear  in  Eq.  (7)  in  terms  of  v±,  E±,  and 
B±.  Moreover,  dfo  {px,Py,Pz,t)  dfo  {Px,Py,Pz,t)  /dpy,  and  dfo  {Px,Py,Pz,t)  jdpx 

must  be  expressed  in  terms  of  5/o  {p±,Pz,C)  /dp±,  dfo  ip±,Pz,0  I^Pz,  dfo  (p±,Pz,|)  /d^, 
because  the  latter  partial  derivatives  are  constant  on  a  characteristic.  Once  these  constant 
derivatives  are  removed  from  the  integral  sign,  the  subsequent  calculation  of  the  Fourier 
transform  of  /j  (z,p,t)  is  facilitated  by  rewriting  them  in  terms  of  dfo  {p±,Pz,<f>,t)  ldp±, 
dfo{p±,Pz,<f>,t)/dpz,  dfo{p±,Pz,<i),t)  /d<f>,  and  dfoipx^PzAi't)  /dt.  The  expression  ob¬ 
tained  for  /i  (z,p,  t)  is 


,  e  PiOc  d  ^ 

+  o  “ITTT'Ht/o 

2  d<j> 


drr  exp  (  i — r  1  Ei-  {z  -f  £  -|-  r) 


U  {p± ,Pz,<f>,t)  J  dr  exp  Ei_  (z  +  £  +  r) 

+  I  exp  {i<f>)  ip±,Pz,  <f>,  t)  J  drr  exp  ^ 

+  V{p±,pz,4>,t)  L  dr  exp  Bi_  (z  -f-  £  -1-  r) 

+  U*  (pj. ,Pz,<i>,t)  J  dr  exp  -^1+  {z  +  VzT,t  +  r ) 

+  \  exp  (-i0)  (PJ- J  exp  f  J  ^i+  +  “^zT, t  +  r) 

+  V*  {p±,pz,<f>,t)  j  drexp^-i^T^  Bi+{z  +  VzT,t  +  T) 


4-e-5— /o 


•/  w  vr- , 

OPz 

epz^c  d 
d(j. 


{Px,Pz,M  [  dTEi^{z  +  v^T,t  +  r) 

J  — OO 

d  f° 

■^fo{p±,Pz,<l>,t)  J  dTTEiz{z  +  VzT,t  +  r) 


In  the  above  equation, 


U  (p±,Pz,  0,*)  =  2  exp 
Tr/_  _  eexp(i<^) 


d  .  d  Pz  d  \ ,  ,  ,  . 

_.p,_  +  )/o  (PX,P.,  <«,i)  . 


In  order  to  relate  components  of  the  perturbed  electromagnetic  fields,  it  is  necessary 
to  obtain  Fourier  (or  Laplace)  transforms  of  fi  (z,p,£)  (Eq.  (17))  and  of  Maxwell’s  equa¬ 
tions  (8)-(ll).  The  simple  assumption  that  all  variables  vary  as  exp  [i  (fcz  —  a;£)]  leads  to 
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inconsistent  results.  In  general,  one  cannot  assume  an  infinite  series  of  terms  of  the  form 
exp  {%  [kz  —  (w  —  nilc/j)  t]}  (where  n  is  an  integer),  because  7  is  not  defined  if  fo  {p±,Pz,0 
contains  an  energy  spread.  In  the  present  analysis,  we  define  the  spatial  and  temporal 
Fourier  transforms  by 

1  poo  poo 

F{k,u))  —  —  I  dz  dtF{z,t)exp[i{u}t  — kz)]^ 

27r  J J—00 

<1  poo  poo 

F{z^t)  =  —  dk  dwF  {k, u)) exp  [i  {kz  —  ujt)] .  (18) 

Stt  J 4/— 00 

The  calculation  of  the  Fourier  transform  of  /i  (2:,p,t),  denoted  by  /i  (A:,p,a)),  requires  use 
of  the  convolution  theorem 

(F1F2)  (Jfc,  w)  =  ^  r  dk'  r  dw'Fi  {k',  u;')  Fi  {k  -k',u-  iv') .  (19) 

27r  J^oQ  J—00 

Then,  using  Eqs.  (17)-(19),  together  with  the  Fourier  transform  of  the  Maxwell  equation 
(9),  and  assuming  that  orders  of  multiple  integrals  can  be  interchanged,  we  obtain  (after 
a  lengthy  calcxilation)  the  following  expression  for  the  Fourier  transform: 

fi{k,p,u)  = 


2  \/2¥ 


<{i(l>)  J 


+00  [  r  n 

du)' <  \u  —  Lo'  —  kvz - - 


'u)-u)'  -kvz  ( .  d  1  5  \  .  Pj.  ck  5  ]  ^  ^ 

V^px  px  90  j a;  -  a;' ^ 


1  e  ,  /•+°°  ^  A  ,  ,  fic 

+--7=  exp (-10)  /  duj  <  u-u  -kVz-i - 

2v27r  J-00  L  7 


’u)-u}' -kv;,  f .  d  ^  I  d  \  ^  .  p±  ck  9  1  ^  ^ 

e  /•+“  f  r  1“^  9 

y  dw'<  u-u’-kv^  i-^foip±,Pz,<f>,<^') 


+  u}  —  u)'  —  kv 


A  \Eu{k,u-J). 
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The  Fourier  transforms  of  the  Maxwell  equations  (8)  and  (10)  are 

/•oo  /^7,‘tc  2 

—  (?k^)  Ei±{k,iij)  =  Amenou)  /  dp±_  /  dpz  / 

Jo  7-00  7o  7^ 

X  exp  (±f<^)  /i  ik,p±,pz,  ,  (21) 

poo  poo  p2n 

kEu{k,u})  =r4nieno  dpj_  dpz  d<f>P±fi  ik,p±,Pz,  <l),u:) .  {22) 

Jo  */— oo  Jo 

By  referring  to  Eq.  (20)  ,  it  is  seen  that  Eqs.  (21)  and  (22)  are  a  set  of  integral  equations 
(with  variable  of  integration  tu')  relating  the  Fourier  transforms  of  the  fields  Ei±  and  Ei^. 

Partial  derivatives  of  fo{p±,Pz,<l>,oj')  appear  in  the  expression  for  fi{k,p,oj)  in 
Eq.  (20).  These  derivatives  are  removed  from  the  integrands  in  Eqs.  (21)  and  (22)  by  inte¬ 
grating  by  parts  with  respect  topj_,  p^,  and  0,  employing  the  periodicity  of  /o  {pXiPz,<l>,(^') 
in  <p.  The  procedure  is  straight  forward  but  requires  much  algebra. 

The  periodicity  of  /o  {p±,Pzi^)  in  ^  (or  <i>)  has  not  yet  been  fully  employed  in  this 
analysis.  Expansion  of  /o  (p±,P2,^)  in  a  Fourier  series  gives 

^  4*cx> 

/o  {p±,Pz,0  =  (Px,Pz)  exp  (tn^) ,  (23) 

^  n=:“00 

where 


1 

Pn  (p±,Pz)  =  c^/o(p±,Pz,^)exp(-m^).  (24) 

Because  /o  {pi,Pz,0  is  real,  g*  {p±,pz)  =  g-n  (pi,Pz).  From  Eq.  (23),  we  see  that  the 
temporal  Fourier  transform  of  fo  {p±,pz,  ^)  =  fo  {p±,Pz,  <f>,t)  is  given  by 


fo  (PJLjPzi^j^) 


+  CX) 

5^  ffn  (p±,Pz)exp(in^)5 
n=— oo 


(25) 


With  the  aid  of  Eq.  (25),  the  integrations  over  <jj'  can  be  carried  out  in  Eqs.  (21)  and 
(22).  Moreover,  with  the  aid  of  the  relation  d<i> exp  {in(f>)  =  27r<5no,  the  integrals  over 
0  can  also  be  completed.  In  fact,  only  the  n  =  0, 1,  and  -1  terms  of  the  infinite  series  in 
Eq.  (25)  contribute  to  the  right-hand  sides  of  Eqs.  (21)  and  (22).  In  the  nonrelativistic 
limit  of 'y  =  1,  the  integral  equations  in  (21)  and  (22)  reduce  to  algebraic  equations  relating 
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the  Fourier  transforms  of  the  fields  Ei±  and  Eiz  because  of  the  Dirac  delta  function  in 
Eq.  (25).  However,  in  the  general  case,  7  =  =  [(p^  +Pz)  * 

Consequently,  Eqs.  (21)  and  (22)  remain  integral  equations  in  the  two  variables  pj.  and  p^, 
relating  such  unknown  functions  as  E\-  (k,iL>  +  IClc/'i  ipxtPz)))  ^iz  +  ^c/l  (PJL>Pi))> 
and  (fc,a;). 

Eqs.  (21)  and  (22)  can  be  reduced  to  a  set  of  integral  equations  in  the  single  vari¬ 
able  p  by  replacing  the  variables  pj.  and  p^  with  the  new  variables  p  and  a,  where 
a  =  tan“^  (j>±/Pz)  is  the  pitch  angle  shown  in  Fig.  1.  Setting  dpxdpz  =  pdpda  in  the 
equations  obtained  from  Eqs.  (21)  and  (22),  we  obtain  the  following  set  of  simultaneous 
integral  equations  relating  the  Fourier  transforms  of  the  perturbed  fields  Ei±  and  Eiz : 


r°°  /  2f2  \ 

D—  {k,u},  Qc)  El-  {k,u)  =  J  dpp^x-+  QcP)  Ei+  ( /!:,a;  -  j 

r°°  f  Q, 

+  dpp\-z{K^,^c,p)Eiz\k,u)- 


f°°  /  2f2c  \ 

D++  {k,  u,  Clc)  Ei+  {k, u))  =  dpp^X+-  {K ^cP)  Ei-  \^k,  uf  -j-  j 

+  dpp^x+z  {k,  a;,  Qc,  p)  Eu  (^k,  u  -f-  ,  (27) 

Dzz  {k,u)  Eu  {k,(j)  =  dpp^Xz+  ik,u),  n^p)  Ei+  (k,uj  - 

+  /  dpp^Xz-{k,u},^lc,p)Ei- (k,u}  + (28) 


In  the  above  integral  equations: 


D —  {k,u,  f2c)  =u}^  -  (?k^  -  J  dp  j  dapgo  (p,  a) 

r2psina  /  kpcosa\  /  fcpcosa 
L  'Y  V  'vm  J  \  “7771 


r2psin<:Y  /  kpcosa\  /  fcpcosa  Qc 
L  7  \  7m  J  \  7m  7 

p^sin^a,  o  o,o\  f  fcpcosa  ^c\~^ 
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Z)-j— I-  fic)  —  ^ ^c)  j 

Dzz  ik,(j)  =  1  -  y/^nl  f  dp  [  dap^go  (p,  a) 

Jo  Jo  I 

(p^cos^a\  f  kpcosa\~^ 

7m  )  ’ 


^  y  2  /  \  sin^a  /  2^0  \ 

X_+ (fcjWjOoP)  =  Jq  \  7  / 


-1 


V  7  /V  7^  7  / 

XH (^>^5  ^C>p)  ~  X— +  (^  5  ^C?p)  > 

X-z{k,u},nc,p)  =  -V^0}nl  [  dapgi{p,a)^;^ 

/  wpcoso;  /  fcp  cosa  Oc\  ^ 

\  7mc  /  V  7m  7  /  ’ 

X+Z  (^>^*^5  ^CjP)  —  X— Z  (^  ■>~^ClP)  1 

'  ,  ^  V  /"  j  ,  vSin^a 

Xz-{k,u;,ac,p) - ^nj(^a;  +  — j  dap^-i  (p,a) 

X  [Pf2!£  L  +  fk')  _  J  L  - 

[  7mc  \  7  /  J  V  7m 

=  5(“+7)'^-(*^'.‘"'  +  7.«-p). 

Xz-f  “  Xz—  (^  1  ^oP)* 


-2 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


(36) 

(37) 


Here,  fip  =  (47rnoe^/m)^/^  is  the  nonrelativistic  plasma  frequency. 

Notice  that  the  structure  of  Eqs.  (26)- (28)  is  a  coupling  of  Ei-{k,u}'), 
Ez  (k,(v'  —  ,  and  Ei+  (k,o}'  —  over  the  range  of  p  for  which  the  equilibrium 

distribution  is  nonvanishing. 


C.  Perturbation  Analysis  for  the  Axial-Dependent  Equilibrium  Distribution. 

The  analysis  for  the  case  of  the  equilibrium  distribution  in  Eq.  (4),  i.e.,  fo  (2,p)  = 
/o  {p±iPz,0  —  fo  (PXjPz>0  “  mQcz/pz)  is  similar  to  the  analysis  of  the  previous  section. 
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The  distribution  is  assumed  periodic  in  C  and  (for  fixed  \z)  in  (j)  with  period  27r.  Conse¬ 
quently,  the  normalization  of  /o  (p±,Pzj  <!>  —  over  momentum  space  is  independent 

of  z.  It  is  defined  to  be  ' 


r  dpz  r 

J —00  J  0 


dp±  d<f>p±fo  \pxypz,<l>-  ==  1- 


The  linearized  Vlasov  equation  for  the  system  is  the  same  as  Eq.  (7)  except  that 
the  faetor  /o  {px)PyiPz,t)  on  the  right-hand  side  is  to  be  replaced  with  fo{z,Px,Py,Pz)- 
Maxwell’s  equations  (8)- (16)  are  applicable  without  modification.  The  derivation  of  in¬ 
tegral  equations  relating  the  Fourier  transforms  of  the  perturbed  fields  involves  a  great 
deal  of  algebra  but  closely  parallels  that  given  in  Sec.  II.B  for  the  spatially-homogenovis 
equilibrium  distribution.  Consequently,  we  omit  the  details  of  the  derivation.  We  remark 
that  in  the  present  derivation  z  plays  much  the  same  role  as  t  in  the  previous  derivation 
and  k  much  the  same  role  as  u).  Periodicity  of  fo  (pj.jP*»C)  hi  C  gives  rise  to  the  Fourier 


senes  expansion 


where 


1 

fo  (p±, Pz,0  =  K  (p±, Pz)  exp  (inC) , 


I 

^nipj-iPz}  J  (P±>Pz>  C)  exp  (  471^).  (■40) 

For  the  z-dependent  equilibrium  distribution  fo{p±,Pzr<P  —  “n^^cz/pz),  the  integral 


equations  are 


2mClt 


DL_{k,u>,Qc)Ei-{k,u)  =  J  dp;,r)-+{k,iJ,Qc,Pz)Ei+ 

+  dp^Ti-z  {k,u,nc,Pz)Bu  (k  +  ^,ui^  ,  (41) 

{k,u},Clc)Ei+  (k,u)  =  J  dp^Tj+^  {k,cj,Q,c,Pz)  Ei-  ^k  - 

+  dpzV+zik,u>,ac,Pz)Eu(^k-^,u^y  (42) 

E'zz  (^>  Elz  (k,  ^)-J  dpg  T]x-  (k,  u,  ClcPz)  El-, 

+  J  dpxVz+ (k,  u,  Qc,  Pz)  Ei+ .  (43) 
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In  the  above  integral  equations: 


-2, 


D'_-  (k,u;,Qc)  =  - J  dPz  dp±b^  (p±,Pz) 

xf2fa-^')L-^-2iy‘ 

L  \  7m  /  \  7m  7  ) 

7^m^cr  ^  ^  \  7m  7  / 

(^5  ^c)  ^  (fc,  61^5  ^c)  ? 

/OO  f^OO 

dpz  /  dp±P±/io(p±,Pz)7~^ 

-00  Jo 


7 


V  72m2c2 ;  V  7^/  ’ 


y/^TT 

77_+  (fc,U),  flcPz)  =  /  clpj.  ^2  (p±,Pi:)  -3-^ 

/Jo  T  ^ 


'y^m^c^ 


X  r_a;2  4.  c2fc  (^fc  +  fa;  -  ^  - 

L  V  Pz  Ji  \  im  nr  J 

77-1 —  ^  P — I- j““^ciPa:)? 

(fc,a;,  fi^Pz)  =  (p±,pz) 

^0 

X  (— a;p2  +  k'fTTK?)  (u)  —  —  — ')  , 

'  \  jm  nr  J 

T7+2  (fc,a;,  HcPr)  =  -^cPz) , 


-2 


Pi 


„  ^  ^  J  ^  ,  ^  pi 

Vz-  {k,uj,  ^cPz)  =  — TT  /  rfpj./i-i  (Pj.,Pz) 


i_ 

'i^rrup' 


X 


fcp^  ^  c^fc 


7m 


^Z+  {kjiO^^C^P^)  Vz—  5  ^ClPz)  • 


(44) 

(45) 


(46) 


(47) 

(48) 


(49) 

(50) 


(52) 


The  structure  of  Eqs.  (41)-(43)  is  the  coupling  of  Ei-{k,u})  to  Build  +  mVlcfPzt^) 
and  Ei+{k'  +  2mQ.c/pz,(j^)  over  the  range  of  p^  for  which  the  equilibrium  distribution  is 
nonvanishing. 
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D.  Some  Special  Cases. 


In  general  Eqs.  (26)-(28)  [Eqs.  (41)-43)]  couple  the  right-  and  left-hand  circularly  po¬ 
larized  radiation  fields  E±  and  the  longitudinal,  relativistic,  plasma  wave  field  E^.  How¬ 
ever,  if  g±i(p±,p^)  =  gii{p,a)  =  0  [h±i(pj.,p^)  =  0],  then  Eq.  (28)  decouples  from 
Eqs.  (26)  and  (27)  [Eq.  (43)  decouples  from  Eqs.  (41)  and  (42)]  to  yield  the  dispersion  re¬ 
lation  for  longitudinal  plasma  oscillations.  The  radiation  fields  remain  coupled,  as  is  seen 

in  Eqs.  (36)  and  (37)  [Eqs.  (51)  and  (52)],  FVom  Eq.  (24)  [Eq.  (40)],  it  is  seen  that  this  sit¬ 
uation  occurs  whenever  the  Fourier  series  for  /o  [fo  0  - 

contains  neither  cos^  nor  sin  ^components.  - 

If  fl±2  (pj,,Pz)  =  g±2  (p.a)  =  0  (Pj.,Pz)  =  0]  then  all  three  fields  and  Eg)  re¬ 
main  coupled,  however  the  radiation  fields  couple  only  through  the  electrostatic  oscillations 
and  not  directly  with  each  other.  Fhom  Eq.  (24)  [Eq.  (40)),  it  is  seen  that  this  situation  oc- 
cum  whenever  the  Fourier  series  for  /„  (px.Pz.tt  -  ^t)  [/„ 
neither  cos  20-  nor  sin  20-components. 

FinaUy,  if  g±iip±,p^)  =  g±2{p±,p^)  =  0  [h±x{pj_,p^)  =  h^^ipx^Pz)  =  0],  then 
Eqs.  (26)-(28)  [Eqs.  (41)-43)]  decouple  completely  and  reduce  to  the  dispersion  relations 

£►__  (/;,w,f2J  =0,  (53) 

D — (A;,a>,  f2c)  =  0,  (5^) 

Dzz{k,u})=Q.  (55) 

These  dispersion  relations  are  identical  to  those  for  the  case  in  which  the  distribution 
m  0  fr^uniformly  random.  Referring  to  Eqs.  (24)  and  (40),  we  see  that  5o(p±,p.)  = 

:^/o  d^fo  (pj.,Pz,0  and  ho  (p±,Pz)  =  Jq''  dqjo  ipx,pz,ri).  For  the  case  of  a  uni¬ 
formly  random  distribution  in  0  (i.  e.,  fo  (px,P.)),  both  of  these  expressions  reduce  to 
9o{p±,Pz)  =  ho{px,Pz)  =  V^h{p±,Pz).  For  this  case  either  Eqs.  (26)  and  (27)  or 
Eqs.  (41)  and  (42)  reduce  to  the  dispersion  relations  for  the  cyclotron  resonance  maser 
with  random  phase  obtained  by  Chu  and  Hirshfield.^ 

The  aualyais  of  Eqa.  (26)-(28)  ((41)-43)]  as  integral  equations  is  the  subject  of  present 
research  and  results  of  the  analysis  will  be  presented  in  a  subsequent  paper.  However, 
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majQy  important  special  cases  exist  in  wHch  the  equations  reduce  to  algebraic  equations 
from  which  dispersion  relations  can  be  derived.  Some  of  these  cases  will  be  analyzed  in 
the  remainder  of  this  paper. 


III.  Dispersion  Characteristics  for  the  Time-Dependent  Equilibrium  Distribu¬ 
tion  without  Energy  Spread. 

In  the  previous  section,  it  was  noted  that,  for  a  uniformly  random  equilibrium  distri¬ 
bution  in  <f>,  equations  (26)- (28)  decouple  and  reduce  to  the  well  known  dispersion  relations 
in  Eqs.  (53)-(55).  In  this  section  it  is  shown  that  as  long  as  there  is  no  energy  spread  in 
the  equilibrium  distribution  in  Eq.  (1),  the  integral  equations  (26)-(28)  reduce  to  algebraic 
relations  between  the  Fourier  components  of  the  fields  even  when  the  distributions  are  not 
uniform  in  (j>.  The  dispersion  characteristics  are  illustrated  with  numerical  examples  for 
the  time-dependent  equilibrium  distribution  with  no  spread  in  pj.  or 

A.  Analysis 

For  a  beam  with  a  definite  energy  7  (po)  rru?  =  'yorru?  =  {pqC^  +  m^c^)  the  most 
general  distribution  in  Eq.  (1)  is  of  the  form 

/o{p,a,0  =  (56) 

P 

where  ^  —  From  Eq.  (24),  the  coeficients  in  the  Fourier  series  expansion  of  this 

distribution  are 

Pn  (p,  a)  =  (p  -  Po)  Pn  (po,  a) , 

(57) 

1  . 

Pn(po,a)  =  -7=  /  d^h{pQ,oc,i)€X.p{-ini) . 

v27r  Jo 

Substitute  Eq.  (57)  into  Eqs.  (29)-(37)  and  then  substitute  the  results  into  the  integral 
equations  (26)-(28).  After  replacing  a»  with  u)  —  2f2c/7o  in  Eq.  (27)  and  with  uj  —  Clc/'yo 
in  Eq.  (28),  we  obtain  three  homogeneous  algebraic  equations  relating  just  three  field 
components.  Expressed  in  matrix  notation,  these  equations  are 

DE  =  0,  (58) 
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where 


D  = 


and 


D — {k,u,nc,po)  .  -x-+{k,u},nc,Po)  -x-z{k,a},nc,PQ)  \ 

-x+-  (k,uj  -  D++  {k,u}  -  -X+z  {k,uj-  ^,f2c,Po) 

\-iCz-{k,uj-^,0.c,P^  -Xz+{k^<^-^,^c,Po^  Dzz{k,w-^,p^  ^ 

(59) 


E  = 


(  Ei-{k,u)  \ 

Eiz(k,(v-^'^ 


(60) 


/ 


The  quantities  Xij  appearing  in  the  matrix  in  Eq.  (59)  are  obtained  from  the  corresponding 
quantities  in  Eqs.  (32)-(37)  simply  by  replacing  each  (p,  a)  with  Qn  (po,or)  and  setting 
p  =  po  and  7  =  70.  Moreover,  from  Eqs.  (29)-(31),  the  diagonal  matrix  elements  can  be 
expressed  as 

, 


D _ {k,u),Clc,po)  =  u)^  —  — 


2  dQgoipo,a) 

r2sina  /  kpocosa\  /  fcpocosa  ^c\ 
1-  70  V  )\  n/om  70/ 


-1 


V  70^ 


'T/qUI^C^ 


a  _  flcN 
70/ 


-2, 


D++  (fc,6t;,nc,Po)  =  D —  (fc,cj,-J2c,po) , 


(61) 


(62) 


Dzz  {k,w, Po)  =  1  -  [  dago  (Po,  oc) 

Jo  7o 

/,  Pocos^a\  /  fcpocosa\~^ 


(63) 


Recall  that,  fip  =  (47re^no/m)^/^  is  the  nonrelativistic  plasma  frequency. 
The  dispersion  relation  for  this  system  is 


detD  {k,u})  =  0. 


(64) 


From  Eq.  (60),  it  is  seen  that,  for  a  given  value  of  the  wavenumber  fc,  the  frequency 
u  is  that  of  the  right-hand  polarized  wave  Ei-.  The  corresponding  frequencies  of  the 
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left-hand  polarized  wave  £?i+  and  of  the  electrostatic  wave  Eu  are  u  -  2ftc/7o  and  u  - 
respectively.  We  point  out  that,  once  the  assumption  of  one-dimensional  spatial 
dependence  is  made  and  equilibrium  self  fields  are  neglected,  the  dispersion  relation  is 
exact  for  equilibrium  distributions  of  definite  energy. 

Again  notice  why,  in  the  general  case,  the  field  components  are  related  by 
integral  equations  (26)-(28)  instead  of  algebraic  equations.  If  the  distribution 
So  {Pi  Oc,<l>  —  ^ct/j(p))  is  nonvanishing  over  a  continuous  range  of  energies  7(p)mc^,  then 
the  mode  Ei-  {k,u})  will  be  coupled  to  a  continuum  of  modes  Ei+  {k,u  —  2flc/7(p))  and 

B.  Case  of  Definite  pj,  and  p^. 


As  numerical  examples,  we  consider  equilibria  of  definite  pj.  =  po  sin  ao  and  Pz  = 
pocosoo-  Consequently,  the  factor  /o  (po,a,^)  in  Eq.  (56)  is 


So  (po,  a,  Ci  =  (^) , 

Sin  a 


(65) 


where  5>  (^)  is  a  function  of  period  27r.  FVom  Eq.  (5),  the  normalization  condition  on  $  (^) 


IS 


(66) 


Moreover,  from  Eq.  (57), 


9n  iPo,oc)  = 

V  27r  sm  a 


(67) 


where 


=  <^^(0  exp  (-171^).  (68) 

Prom  Eq.  (66),  sq  =  1.  Also  notice  that  =  s*. 

Substituting  Eq..  (67)  into  Eqs.  (61)-(63)  and  into  Eqs.  (29)-(31)  with  (p,  a)  replaced 
by  ffn  (Po)  oc) ,  we  can  express  the  dispersion  relation  in  Eq.  (64)  as  the  following  tenth  degree 
polynomial  equation  (in  either  Cj  or  k)  with  real  coefficients: 

M__  A4+ = 
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1  f 

2  (^)  (k,u,)  +  -i)  V..  ft, a)] 


(69) 


1  l'4 


4  /»!  (»2»ix + (i,»  -  aD  -  *2)  ^ 

and  k  -  ck/uj^.  Dimensionless  velocities  are  given  bv  5.  =  ,,  /  j  ^  ^ 

=p„sh..o/.(;.)„.  and  =p„cos„„/,(p„,„.  ““ 

(*.4  =  (ch  -  -  l) ^  ^ 

-  1)  ^  ^  (1  - /?f )  . 

/  a;2  ''  (72) 

The  dispersion  relation  in  Eq.  (69),  which  is  valid  for  both  eomnl  -  a 

is  invariant  under  the  transfomration  "  “ 

k  — ^  -k\ 

U)  — >  -cb*  +  2. 

(73) 

Tv  Tr“"  ‘  ““’  '•""“  -  ■»'’“■  “• —I"-.- » T  r 

ments  of  D  and  E  in  Eq.  (58).  If  either  sr  ^  0  or  s,  ^  0,  then 

M^:i:f)  .  E^(b,rh) 

^‘-(*’")  ^l-  (-**.-<h*+2)  Er^.(i,a,_2)’ 
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and,  if  5i  0, 


Ei+(k,u-2)  Ei+(-k*,-u;*)  El.(k,u) 

Eu  {k  w  - 1)  Eu  {-y,  -1^*  + 1)  El,  (k,  d)  -  l)  ’ 

El-  (k,  d;)  El-  (-k*,  -u*  +  2)  _  El^(k,u-2) 

Eu  (k,u>-l)  Eu  {-k*,  -w*  + 1)  El,  (k,u>-l)' 


(75) 


Equation  (69)  gives  the  ten  branches  of  the  dispersion  relation  u){k).  -The  behavior  of 
C}/h4.0T  large  \k\  is  easily  determined  for  each  of  these  branches.  As  |^|  — >  oo,  u>{k)/k  — >  +1 
for  two  branches,  C}{k')/k  — >  — 1  for  two  branches,  and  Cb(k)/k  — *  Pz  for  six  branches. 

Simple  expressions  axe  easily  obteiined  for  the  large  |fc|  behaviors  of  d>(fc)  for  all  ten 
branches  if  either  si  or  S2  vanishes  in  Eq.  (68).  If  Si  =  0  and  S2  ^  0,  then  [from  Eqs.  (57) 
and  (67)]  gi  (j>,Oi)  =0  and  p2  {p,ct)  0-  If  follows  from  the  discussion  in  Sec.  II. D  that 
the  electromagnetic  components  are  coupled  and  the  electrostatic  component  is  imcoupled. 
[Such  a  situation  holds  for  (but  is  not  exclusive  to)  the  distribution  $(0  in  Fig.  3(b), 
provided  that  the  parameter  a  ^  tt.  For  this  distribution,  Si  =  0  and  S2  =  sin  a/a. ]  Two 
of  the  branches  pertain  to  the  uncoupled  electrostatic  waves  and  obey  the  exact  dispersion 
relations 

u  =  kl3z  +  l±'^0--Pl)^^^  ■  (76) 

UJc 

for  sufficiently  large  |fc|,  the  remaining  eight  branches  obey  the  approximate  dispersion 
relations 


d>  Ci  icPz  +  1  ±  (1  +  lS2l)^^^  . 

v2 

(77) 

u)^kl3z  +  l±-^—p±(l-\s2\)^^^, 

(78) 

/  2\1/2 
+  , 

(79) 

(80) 
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For  real  k,  Eqs.  (77)  and  (78)  give  two  branches  with  positive  Ima>  provided  that  jsal  <  1- 
In  this  case,  growth-rate  curves  (Imw  vs.  real  k)  will  show  two  unstable  branches  at  large 


k. 

Usi^O  and  S2  =  0,  then  [from  Eqs.  (57)  and  (67)]  (p,  a)  =  0  and  gi  (p,  a)  #  0.  It 

follows  from  the  discussion  in  Sec.  II.D  that  the  transverse  electromagnetic  waves  are  cou¬ 
pled  through  the  longitudinal  electrostatic  wave.  [A  nonexclusive  example  of  a  distribution 
$(^)  having  this  property  (if  K"  7^  0)  is  presented  in  Fig.  3(c).  Using  Eq.  (68),  we  find  that 
—  -4iK  and  S2  =  0  for  this  example.]  In  this  case,  the  large  Ifcj  approximations  for  four 
of  the  ten  branches  of  the  dispersion  relation  are  the  same  as  those  given  m  Eqs.  (79)  and 
(80)  for  the  distribution  previous  case.  Approximations  for  the  remaining  sbc  branches  are 


[The  maximum  possible  value  of  AT  in  Fig.  3(c)  is  l/27r,  and  the  corresponding  maximum 
value  of  |sil  is  2/tt  =  0.6366....  It  is  evident  from  Eq.  (84)  that  the  k±  are  real  and  positive 
for  all  isil  <  I/V2  =  0.7071....  Consequently,  Eqs.  (81)  and  (83)  provide  for  two  unstable 
modes  at  large  values  of  real  k  for  the  distribution  in  Fig  4(c).] 

A  nonexclusive  example  of  a  distribution  for  which  neither  Si  nor  S2  vanishes  (unless 
the  parameter  a  =  tt  or  27r)  is  presented  in  Fig.  3(d).  From  Eq.  (68),  it  follows  that  in  this 
example  si  =  (2/a)  exp  (-ia/2)  sin  (a/2)  and  S2  =  (1/a)  exp  (-ia)  sin  a.  If  both  5i  and  S2 
are  nonvanishing,  determining  the  large  |fc|  behavior  of  the  dispersion  relation  in  Eq.  (69) 
is  more  difficult  than  in  the  previous  cases.  Four  of  the  large-|fc|  branches  are  given  by 
Eqs.  (79)  and  (80).  The  behaviors  of  the  remaining  six  branches  (including  all  that  may 
show  growth  at  large,  real  k)  are  determined  by  solving  a  cubic  equation  in  {u)-kl3^  -  1)^. 
Further  details  will  not  be  given  in  this  paper. 
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Unless  Si  =  S2  ==  0,  some  of  the  eigeamodes  E  ia  Eq.  (58)  will  involve  two  or  more  of 
the  components  E-{k,u)),  E+{k,u-2),  and  Eg{k,C)- 1).  A  quantity  that  will  be  employed 
to  measure  the  relative  importance  of  E-{k,Co)  and  E+{k,u)  —  2)  is  the  Poynting  flux  ratio 
defined  by 

This  quantity  is  the  ratio  of  the  time-averaged  Poynting  vectors  that  the  field  associated 
with  each  component  would  produce  in  the  absence  of  the  other  component.  A  time- 
dependent  interference  term  due  to  the  different  frequencies  of  the  components  is  not 
included. 

C.  Numerical  Examples. 

In  the  following  numerical  examples,  k  is  restricted  to  be  real.  Then,  ImtD  >  0  indicates 
an  unstable  mode.  If  ^  is  restricted  to  be  real,  then  the  transformation  in  Eq.  (73)  is 
equivalent  to  inverting  a  plot  of  RecD  vs.  k  (real)  through  the  point  (k,Reu^  =  (0, 1)  and 
reflecting  a  plot  of  ImcJ  vs.  k  (real)  through  the  Ima>-axis.  It  follows  from  the  invariance  of 
the  dispersion  relation  in  Eq.  (69)  under  this  transformation  and  from  Eqs.  (74)-(75)  that 
there  is  no  loss  of  generality  if  numerical  examples  are  limited  to  the  case  of  nonnegative 
real  k. 

Parameter  values  in  all  of  the  numerical  examples  below  in  Figs.  4-7  are  Oplu)^  =  0.05, 
7o  =  2,  and  Qq  =  0.4.  In  order  to  ensure  that  values  selected  for  Si  and  S2  are  realistic 
(i.e.,  correspond  to  $(^)  >  0  in  Eq.  (65)),  we  eissume  that  4>(^)  has  one  of  the  functional 
forms  shown  in  Figs.  3(a)-3(d). 

Example  1:  If  the  distribution  $(^)  is  uniform  [see  Fig.  3(a)],  then  si  =  S2  =  0. 
[Such  a  distribution  is  also  attained  with  a  =  tt  in  Fig.  3(b),  K  =  0  in  Fig.  3(c),  or 
o  =  27r  in  Fig.  3(d).]  In  this  case,  the  dispersion  relation  in  Eq.  (69)  decouples  into 
the  three  independent  relations  M —  (k,u}^  =  0,  M++  (k,u^  =  0,  and  Mgz  = 

0.  These  are  respectively  the  dispersion  relations  for  uncoupled  right-  arid  left-handed 
circularly  polarized  electromagnetic  waves  and  for  the  electrostatic  wave.  For  a  given  k, 
the  frequencies  of  these  waves  are  u,  u  —  2,  and  th  —  1,  respectively.  Growth-rate  curves 
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Ac^»ntpfl  in  Fic.  4{a)  over  the  interval  0  <  fc  ^ 

(Imw  vs.  k)  for  this  limiting  case  are  pr  ^  ^  presented  in  Fig.  4(b). 

Corresponding  plots  of  R«h  vs.  1=  over  the  m  -  ^  ^ 

Letters  on  these  plots  designate  correspon  obtained 

d/  X  •va.irIP  The  erowtb-rate  curve  segments  CL»G  ana  ni 

and  C  in  Fig.  4(a)  comcide.  Ibe  gr  ,  «  .  (r^R)  has  only  E.{k,u)  as 

r  M  ffc  u;"!  =  0,  and  a  corresponding  eigenmode  B  m  Eq.  (58) 

from  M-  J»=.u-)  .  are  obtained 

a  nonvanishing  ^mponent.  gr  ^  pppvanishing  component, 

from  Af++ (b.d-)  =  0- A  conrespondmg  e, genmode  has 

namely  -  2).  _  ^  „tation  of  the 

The  plot  of  Rew  vs.  k  in  Fig.  U  „  /f  p.Jik  u)  are 

,n.h.rate  corves  in  Fig.  4(a,.  ff^  >  0  “;:^;^^tvesl:  travel  in  the  iorwmd 

(positive-z)  doectron.  i  ^  gbnilarily, 

is  unchanged,  however  they  if  fc  >  0  and  Bed,  <  0. 

the  electrostatic  wave  associated  wi  z(  i  )  4 r  ^  ronsider  the  segmented 

4-  f  th.  use  of  Fig  4(b)  in  interpreting  Fig.  4(a),  consider  the  segm 
As  an  illustration  of  the  use  ot  g.  i  ;  .  r  cvclotron  maser  insta- 

A(nr\r)r  which  gives  the  growth  rate  of  the  cyclotron  m 
g.owth-rate  curve  A{BC)DG,  S  ^  ^  Ref^  _  2  <  0 

bility.  Segment  AB  pertains  to  +{ 

j  X  i;«rr  T  WP  electron! oimil  y>  & 

.„g  backwar  -  rave  ^  ^  ^  ^  ^  ^  evorywhcre  on  CDG  so 

pertarns  to  E.  (k,  u,).  forward-travehng, 

,bat  growth-ra.  cmve  ^  ^  ,be  growth-rate  curve  ff  f  for  the 

RHP  clectromagnetrc  waves.  Srm  V  electromagnetic  waves.  Also,  the 

whistler  instabUity  pertains  to  forward-trave  mg,  ^ 

4  17  F  nprtains  to  backward-traveling  LHP  waves,  ana  b 

grewth-rate  curve  segment  ^  ^be  roots  of  (fc,-) 

reentEJ  pertains  to  forward-travelmg  LHP  waves.  tJ 
are  real  no  corresponding  growth-rate  curves  appear  m  Fig.  (a). 

Electromagnetic  and  beam  waves  axe  said  to  be  in  ain  irgiven  by 

instability  when  d,  =  k  and  d,  =  fcfe  +  L  the  well-known 

\  ■'r  -  f-  -  1/(1  -  Pz)-  In  tills  example,  kr  -  4.94.  rig.  ‘iv  ; 

\  oJr-kr-i-IK  ^  ^kU-k  in  an  uncoupled  system. 

\  fact  that  no  growth  of  RHP  radiation  occurs  at  fc  -  k. 
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Finally,  it  is  emphasized  that  no  special  relation  exists  between  the  LHP  and  RHP 
waves  considered  above  when  si  ■=  sj  However,  the  plots  in  Fig.  4  will  be  approached 
by  any  system  using  our  parameters  in  the  limit  in  which  both  Si  and  S2  approach  zero. 

Example  2:  An  explicit  example  of  a  nonuniform  distribution  in  phase  is  obtained 
by  selecting  o  =  7r/4  in  Fig.  3(b).  Then  the  unperturbed  electron  beam  consists  of  two 
streams  with  respective  distributions  centered  about  ^  =  0  and  ^  =  tt.  Each  distribution 
is  a  water  bag  of  width  7r/4.  The  corresponding  parameters  defined  in  Eq.  (68)  are  Si  =  0 
and  S2  =  2\/^/7r,  Because  Si  =  0,  the  eigenmodes  E  in  Eq.  (58)  are  of  two  types.  The  first 
type  of  eigenmode  has  only  Eiz{k,u))  as  a  nonveinishing  component.  The  corresponding 
dispersion  relation  is  Mzz{k,u)  =  0,  which  does  not  allow  for  growth.  [See  Eq.  (76).] 
The  second  type  of  eigenmode  has  two  nonvanishing  components,  namely  Ei-{k,u)  and 

£^14- (fc, —  2).  Some  of  these  eigenmodes  are  unstable. 

Growth-rate  curves  for  this  system  for  0  <  ^  <  15  and  corresponding  plots  of  RetD 
Vi.  k  (for  0  <k  <  1.5)  are  presented  in  Figs.  5(a)  and  5(b).  As  a  measure  of  the  relative 
importance  of  the  RHS  and  LHS  electromagnetic  waves  associated  with  unstable  modes, 
plots  of  the  Poynting  flux  ratio  in  Eq.  (85)  as  a  function  of  k  (for  unstable  modes  only)  are 
presented  in  Fig.  5(c).  Letters  show  corresponding  points  on  Figs.  5(a)-5(c).  By  comparing 
Figs.  5(a)  and  5(b),  we  see  that  the  growth-rate  curve  segments  BA,  CD,  and  FE  pertain 
to  xxiodes  consisting  of  a  forward-traveling  RHP  electromagnetic  wave  (because  Rew  >  0) 
and  a  backward-traveling  LHP  electromagnetic  wave  (because  Rea)  —  2  <  0).  All  other 
segments  of  the  growth-rate  curves  pertain  to  modes  consisting  of  forward  traveling  RHP 
and  LHP  electromagnetic  waves. 

In  the  case  of  a  uniform  distribution  in  ^  (Example  1),  no  growth  of  RHP  electro¬ 
magnetic  waves  occurs  at  the  resonance  wavenumber  A:  =  fcr  =  4.94.  [See  Fig.  4(a).] 
The  growth-rate  curve  CDH  in  Fig.  5(a)  shows  a  mode  a.t  k  —  kr  ~  4.94  which  grows 
significantly  faster  than  any  mode  in  Fig.  4(a).  From  Fig.  5(a),  we  see  that  the  RHP 
Poynting  flux  associated  with  the  mode  is  almost  twenty  times  the  LHP  Poynting  flux. 
Consequently,  growth  of  RHP  radiation  is  now  possible  at  k  =  kr  although  if  must  be 
accompainied  by  a  smaller  growing  component  of  LHP  radiation.  We  remark  that  among 
computations  so  far  carried  out  those  for  systems  with  two-stream  distributions  in  ^  (with 
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a  phase  difference  of  tt)  show  the  most  rapid  growth  rates  a.t  k  =  kr. 

Referring  to  Fig.  5(a),  we  see  that  there  are  two  unstable  branches  at  large  k  in 
conformity  with  Eqs.  (77)-(78).  From  Fig.  5(c),  it  is  seen  that,  for  either  branch  at  large 
k,  the  RHP  Poynting  flux  is  approximately  double  that  of  the  LHP  Poynting  flux.  Finally, 
note  from  Fig.  5(c)  that  for  unstable  eigenmodes  at  small  k  the  ratio  of  the  backward 
traveling  LHP  Poynting  fliot  to  the  forward  traveling  RHP  Poynting  flux  depends  very 
strongly  on  the  branch  of  the  dispersion  relation  and  varies  rapidly  with  k  for  a  given 
branch. 

Example  3:  As  a  second  example  of  a  nonuniform  phase  distribution,  select  the  form 
of  $(^)  in  Fig.  3(c)  and  choose  the  parameter  value  K  =  l/27r.  [Equivalently,  we  could  let 
a  =  TT  in  Fig.  3(d).]  Then  the  equilibrium  particle  phases  are  uniformly  distributed  between 
^  =  0  and  ^  =  TT,  and  no  particles  have  phases  in  the  range  tt  <  |  <  27r.  Fourier  components 
in  Eq.  (68)  are  Si  =  —2ifTr  and  S2  =  0.  Because  $2  =  0,  the  field  components  Ei-{k,u)), 
and  Ei+{k,u)—2)  are  indirectly  coupled  through  the  electrostatic  component  Eiz{k,i^  —  1)- 
Consequently,  the  eigenmodes  E  in  Eq.  (58)  will  (in  general)  have  three  nonvanishing  field 
components.  Growth-rate  curves  for  this  example  are  presented  in  Fig.  6(a)  for  0  <  fc  <  10, 
and  corresponding  plots  of  Rea)  vs.  k  (for  the  interval  0  <  <  2)  are  given  in  Fig.  6(b). 

Plots  of  the  Poynting  flux  ratio  in  Eq.  (85)  are  shown  in  Fig.  6(c).  Finally,  as  a  measure 
of  the  relative  importance  of  the  electrostatic  component  of  the  unstable  eigenmodes,  we 
present  a  plot  of  2~^^^\Ei-{k,u})/Eiz{k,u)  —  1)|  vs.  k  in  Fig.  6(d).  (The  factor  of  2~^/^ 
appears  in  the  field  ratio  because  2~^/^E±  is  the  proper  normalization  of  coefficients  of 
the  complex  unit  vectors  for  LHP  and  RHP  waves  when  comparison  is  to  be  made  with 
Cartesian  field  components.)  Letters  show  corresponding  points  in  Figs.  6(a)-6(d). 

Reference  to  Figs.  6(a)  and  6(b)  shows  that  the  wave  associated  with  the  compo¬ 
nent  Ei+{k,oj  —  2)  is  left-hand  polarized  and  backwards  traveling  for  eigenmodes  on  the 
growth-rate  curve  segments  GH,  ABC,  and  DEF.  The  electrostatic  wave  associated  with 
Eiz(k,u}  —  1)  is  forward  traveling  for  all  unstable  modes  except  for  those  modes  on  the 
growth-rate  curve  ABC  for  which  k  is  very  close  to  zero.  All  other  components  of  unstable 
eigenmodes  represent  forward-traveling  waves. 

Referring  to  Fig.  6(a),  we  see  that  two  unstable  branches  of  the  dispersion  relation 
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are  present  at  large  values  of  in  conformity  with  Eq.  (81)  and  (83).  BVom  Figs.  6(c) 
and  6(d),  it  is  seen  that  the  RHP  Poynting  flux  exceeds  that  of  the  LHP  electromagnetic 
wave  by  a  factor  of  approximately  three  for  both  branches  at  large  k.  The  electrostatic 
contribution  to  eigenmodes  on  the  upper  branch  is  relatively  very  small.  On  the  other 
hand,  the  electrostatic  fleld  amplitude  in  eigenmodes  on  the  lower  branch  is  of  the  same 
order  of  magnitude  as  the  LHP  electromagnetic  fleld  amplitude. 

Two  branches,  MNO  eind  GHI,  show  moderate  growth  rates  at  the  resonance 
k  =  kr  =  4.94.  Reference  to  Figs.  6(c)  and  6(d)  shows  that  the  RHP  Poynting  flux  is 
significantly  larger  than  the  LHP  Poynting  flux  for  the  eigenmodes  associated  with  either 
of  these  branches  at  ^  =  fc,..  Moreover,  \Ei—(k,u)')/Eiz{k,Cij  —  1)|  >  10  for  either  branch 
at  k  =  kj..  Again,  this  behavior  is  in  contreist  with  the  case  of  a  uniform  distribution  in  ^ 
where  no  growth  of  RHP  electromagnetic  radiation  takes  place  at  the  resonance  value  of 
k.  However,  the  growth-rates  at  k  =  kr  in  Fig.  6(a)  are  not  large,  being  slightly  less  than 
the  maximum  growth-rates  that  appear  in  Fig.  4(a)  for  the  case  of  a  uniform  distribution. 

Although  the  growth-rate  peak  ABC  is  very  narrow,  eigenmodes  at  points  near  its 
maximum  are  the  fastest  growing  modes  of  this  system.  Moreover,  reference  to  Fig.  6(c) 
shows  that  these  modes  contain  a  relatively  strong  backward- traveling,  LHP  component. 

Example  4:  As  our  final  numerical  example,  we  treat  the  limit  of  o  =  0  for  the  distri¬ 
bution  in  Fig.  3(d).  In  the  limit,  the  distribution  becomes  $  (^)  =  Y,-oo  2«7r)  with 
Si  =  S2  =  1-  In  this  case,  ^  =  Uct  for  all  particles  in  the  equilibrium  beam.  Growth-rate 
curves  for  the  interval  0  <  ^  <  8  and  plots  of  Rea>  vs.  k  for  the  interval  0  <  <  1.5  appear 

in  Figs.  7(a)  and  7(b),  respectively.  For  unstable  eigenmodes,  plots  of  the  Poynting  flux 
ratio  in  Eq.  (85)  vs.  k  and  2~^^^\E\^{k^Cj) / Eiz{k,u  — 1)|  vs.  k  are  presented  in  Figs.  7(c) 
and  7(d),  respectively.  Letters  on  these  graphs  show  corresponding  points.  Eigenmodes 
belonging  to  growth-rate  curve  segments  ABC  and  DE  have  backward-traveling  LHP 
components  and  forward-traveling  RHP  and  electrostatic  components.  A  tiny  growth-rate 
peak  appears  at  G  in  Fig.  7(a).  Reference  to  Fig.  7(b)  shows  that  its  LHP  and  electro¬ 
static  components  are  backward  traveling.  Eigenmodes  on  all  other  segments  contain  only 
forward-traveling  components. 

A  striking  feature  of  the  growth-rate  curves  in  Fig.  7(a)  is  that  no  growth  occurs 
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for  values  of  k  greater  than  approximately  6.2.  That  is,  no  branches  of  the  dispersion 
relation  show  growth  in  the  limit  of  large  k,  [Suppression  of  instability  in  the  whistler  by 
the  electrostatic  wave  is  disciissed  in  Ref.  9.]  Also,  notice  the  interval  of  no  growth  FH 
(1.4 -fc -2). 

Another  striking  feature  is  the  great  height  of  the  growth  peak  ABC.  Prom  Fig.  7(c), 
we  see  that  the  Poynting  flux  of  the  backward-traveling  LHP  electromagnetic  wave  ex¬ 
ceeds  that  of  the  forward-traveling  RHP  electromagnetic  wave  over  most  of  the  interval 
of  this  growth  peak.  The  amplitude  of  the  forward-traveling  electrostatic  wave  is  seen 
[from  Figs.  7(c)  and  7(d)]  to  be  of  the  same  order  of  magnitude  as  the  amplitudes  of  the 
electromagnetic  waves. 

Again,  we  see  growth  of  RHP  electromagnetic  waves  at  the  resonance  wavenumber 
_  4  94,  In  fact,  the  maximum  of  the  growth-rate  curve  HJJ  in  Fig.  7(a)  is  situated 
very  close  to  the  resonance  wavenumber,  and  the  growth  rate  at  this  maximum  exceeds 
any  growth  rate  for  the  gyrotropic  case  is  Fig.  4(a).  From  Figs.  7(c)  and  7(d),  it  is  seen 
that  the  largest  component  for  eigenmodes  near  this  maximum  is  that  corresponding  to 

forward-traveling  RHP  electromagnetic  radiation. 

To  summarize,  it  is  evident  that  a  richness  of  structure  in  the  growth-rate  curves  can 
be  produced  by  introducing  nonuniform  distributions  $(0-  Using  proper  choices  of  $(^), 
temporal  growth  rates  near  k  =  kr  can  be  significantly  increased  and  growth  rates  at  large  k 
can  be  on  the  one  hand  enhanced  or  on  the  other  hand  completely  suppressed.  Two-stream 
equilibrium  distributions  such  as  that  in  Example  2  seem  to  be  most  effective  in  enhancing 
growth  rates  at  the  resonance  k  =  kr-  Gaps  of  no  growth  can  be  introduced  at  moderate 
values  of  k  At  small  values  of  k  where  eigenmodes  may  contain  backward-traveling  compo¬ 
nents,  growth  rates  and  the  properties  of  eigenmodes  can  be  greatly  changed  by  changing 
$(^).  [This  latter  fact  suggests,  but  does  not  prove,  that  absolute  instability  properties 
may  depend  strongly  on  4>(4).  However,  no  pinch-point  analyses  of  these  systems  have 

been  carried  out.^®"^®] 
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IV.  Dispersion  Characteristics  for  the  Axial-Dependent  Equilibrium  Distribu¬ 

tion  without  Ajcial  Momentum  Spread. 


A.  Analysis 

Finally,  we  consider  the  axial-dependent  equilibrium  distribution  fo{p±,Pz,0  — 
fo  {p±^Pz,<f>  ~  mQ,czlpz)  for  which  the  perturbed  field  components  are  related  by  the  in¬ 
tegral  equations  (41)-(43).  Even  if  the  equilibriiim  distribution  is  not  uniformly  random 
in  <i>,  Eqs.  (41)-(43)  will  reduce  to  algebraic  equations  if  there  is  no  spread  in  in  the 
equilibrium  distribution.  (A"spread  in  energy  is  permitted  if  it  is  due  only  to  a  spread  in 
px*)  The  most  general  equilibrium  distribution  having  this  property  is 

foip±,Pz,0=^iPz-Pzo)foip±,PzO,0-  (86) 

It  follows  from  Eq.  (40)  that 

{p±iPz)  ~  ^  {Pz  Pzo)  hn  (PJ-jPzO)  )  (87) 

where 

1 

Ktip±,Pzo)  =  -^  dC/o(p±,PzO,Oexp(-mC).  (88) 

vStt  Jo 

Three  homogeneous  equations  relating  just  three  field  components  are  found  using  a  pro¬ 
cedure  similar  to  that  used  in  obtmning  Eq.  (58).  In  matrix  form  the  equations  are 

D'E'  =  0,  (89) 


where 


D'  =  - 


DL-ik,u,Qc,Pzo)  -rj-+{k,u,^lc,Pzo). 

-V+-  [k  -t-  ^^,u;,Qc,Pzo)  D'++  (k  +  ^^,a;,fic,Pzo) 

-^z-  (k  +  ^yU^QcPzo)  -fiz+  (k  -I-  ^,a;,f2c,p^o) 


-rj-z  ik,u,Clc,Pzo) 
~V+z  (k  -h  fioPzO 

(90) 


E'  = 


El-  {k,u) 


Elz  {k  + 
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The  quantities  f}ij  in  Eq.  (90)  axe  obtained  from  the  corresponding  quantities  in  Eqs.  (47)- 
(52)  by  replacing  (pj.,Pz)  with  hn  (Pi.,Pz)),  Pz  with  p^o,  and  7  with  7(pj.,Pzo)-  The 
remaining  quantities  in  Eq.  (90)  can  be  obtained  from  Eqs.  (44)- (46).  They  are 


DL^  (fc,  w,  J^cPzo)  =  a;2  -  - dp±ho  (px» Pzo) 

xhfu-  ^  (^J  _  fcPzO  fic 

I-  V  7(PJ.,Pzo)my  V  7(P±>Pzo)^w  7(p±.pzc 

-  ,,  (O.^  -  (■.' - - ,  ,  (92) 

T  {p±,pz0)  ''  V  7  {p±,PzO)  m  7  (p±, P20)  J  -I 


(fc, a;,  f2c)Pzo) — -D — (/c,  Ci^,  f2cjPzo)) 


(93) 


roo 

(^.  ‘^^.Pzo)  =  1  -  J  dpxpxho  (px.Pzo)  7~^  (Px,Pzo) 

o(“-555fcr- 


X  1 


P*o 


'r^(p±,PzQ)m^c^J  'l({p±,PzQ)m^ 

Once  the  assumption  is  made  that  equilibrium  self  fields  can  be  neglected,  the  exact 
dispersion  relation  for  the  case  of  definite  Pz  =  PzO  is 


detD'  (A:,  a;)  =  0. 


(95) 


For  a  given  frequency  u,  the  wavenumber  of  the  right-hand  polarized  wave  Ei-  is  k.  The 
wavenumbers  of  the  fields  Ei+  and  Eiz  are  k  -f-  2mQclPzO  and  k  +  mQc/PzO,  respectively. 
This  result  shows  why  mtegr£d  equations  relate  the  field  components  if  a  spread  in  pz  is 
present  in  fo{p±,Pz,  <f>—mQcz/pz)-  Then,  Ei-  {k,u})  couples  to  Ei+  {k  +  2mQ,c/pz,(^)  and 
to  Eiz  {k  -f  mTlc/Pz^i^)  over  a  continuum  of  values  of  mClc/Pz" 

Finally,  we  emphasize  that  the  eigenmode  E  in  Eq.  (60)  is  of  a  different  nature  than 
the  eigenmode  E'  in  Eq.  (92).  The  eigenmode  E  is  a  composite  of  three  components  which 
refer  to  waves  of  the  same  propagation  vector  but  of  different  frequencies.  These  frequencies 
differ  by  fixed  real  values.  On  the  other  hand,  the  eigenmode  E'  is  a  composite  of  three 
modes  which  refer  to  waves  of  the  same  frequency  but  of  different  ..propagation  vectors. 
These  propagation  vectors  differ  by  fixed  real  values.  In  either  case,  if  is  restricted  to 
real  values,  then  temporal  growth  or  decay  rates  are  given  by  Ima>  for  all  components.  If 
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Cb  is  restricted  to  real  values,  then  spatial  growth  or  decay  rates  are  given  by  Imfc  for  all 
components. 

B.  Case  of  Definite  and  pz 

As  in  Sec.  Ill,  we  present  numerical  examples  for  cases  in  which  both  Pz  =  PzO  and 
P±  =  PXQ  have  definite  values  in  the  equilibrium  distribution.  Consequently,  the  factor 
/o  (p±>PiO>C)  bi  Eq.  (86)  is  selected  to  be  of  the  form 

So  C)  =  (C) ,  -  (96) 

Px 

where  'J'  (C)  is  a  periodic  function  of  C  (or  (j))  of  period  27r.  BVom  Eq.  (5),  the  normalization 
condition  on  'I'  (C)  is 

d(t>^  dC^  (C)  =  1.  (97) 

It  follows  from  Eq.  (88)  that 

hn{px,Pzo)  =  (98) 

V27r  Px 

where 

y^n-  [  (C)  exp  (-inC) .  (99) 

Jo 

Notice  that  =  1  and  that  w-n  =  u;*. 

With  the  aid  of  Eqs.  (92)-(94)  and  Eqs.  (46)-(52),  we  obtain  the  dispersion  relation 
in  Eq.  (95)  for  the  case  of  definite  p±  and  pz.  The  dispersion  relation  is 

M--(k,<:b)M!^.+  (k,u)Mzz[k<^)  = 
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Dimensionless  frequencies,  wavenumbers,  eind  velocity  components  are  defined  as  in 
Sec.  III.B,  and 


4-  J  -  l)  -  ^  -  2)  1) 

(u|)].  (101) 


The  quantities  M—  ^k,u>J  and  Mzz  ^k,u>^  are  defined  in  Eqs.  (70)  and  (72),  respectively. 

T^jWp.  Eq.  (69),  Eq.  (100)  is  valid  for  complex  ui  and  complex  fc.  It  is  invariant  under 
the  transformation 


— >  —(2)*. 


(102) 


In  analogy  with  Eq.  (74),  it  follows  from  the  matrix  equation  (90)  that,  if  either  Wi  0 
or  W2  7^  0,  then  under  the  transformation  in  Eq.  (102) 

Ei+ +  2//3z,‘^^  E*_^k,u)^  (103) 

El-  (fc,  (2i)  El-  (^-k*  -  2l^z.  -w*)  Et+  (k  +  2//3„  ni)  ’ 

Moreover,  in  analogy  with  Eq.  (75),  if  ini  71^  0,  then  under  the  transformation  in  Eq.  (102) 

Ei+  {k  +  2/Pz,^y  Ei+{-k\-u*^  El-  (k,Q) 

Eizik  +  l/Pz^C;)  Eiz{-k*-l/fSz,-C^*)  EJ,  (fc  +  l/)0„ w)  ’ 

(104) 

El-  (k,u)  El-  (-k*  -  2/ IS,,  -to*)  ^  EI+  (k  +  2/l3z,tb) 

Eiz  (k  +  l/l3z,Co)  Eiz  {-y  -  l/lSz,  -u^*)  El,  (k  +  l/l3„Cb)  ' 

The  ratio  of  the  time-averaged  Poynting  vectors  associated  with  the  individual 
Ei+  ^k  +  2/fiz,cb)  and  Ei-  fields  is 
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Like  the  dispersion  relation  in  Eq.  (69)  of  Sec.  III.B,  Eq.  (100)  is  a  tenth-degree 
polynomial  equation  in  either  k  or  Cb.  The  large  |fc|  behavior  of  u){k)/k  given  by  Eq.  (100) 
is  the  same  as  that  given  by  Eq.  (69),  including  two  branches  with  C}{k)/k  a  1,  two 
branches  with  u>{k)/k  Ci  —1,  and  six  branches  with  u){k)/k  ci  kfiz- 

If  either  wi  or  vanishes,  approximations  for  C){k)  valid  for  large  values  of  k  are 
readily  determined  and  axe  found  to  be  very  similar  to  those  found  in  the  previous  section 
for  the  case  of  the  time-dependent  equilibrium.  If  toi  =  0  and  W2  ^  0,  then  the  transverse 
electromagnetic  components  are  coupled  to  each  other  and  the  electrostatic  component 
is  uncoupled.  The  large  |^|  behaviors  of  the  ten  branches  of  the  dispersion  relation  in 
Eq.  (100)  are  similar  to  those  given  in  Eqs.  (76)- (80)  for  the  corresponding  case  (si  = 
0,  S2  7^  0)  of  the  dispersion  relation  in  Eq.  (69).  Four  of  the  branches  now  obey  Eqs.  (76) 
£ind  (79).  Four  branches  obey  Eqs.  (77)  and  (78)  with  §2  replaced  by  W2,  that  is 

u^kl3,  +  l±±‘^p^{l  +  \w2\)'/\  (106) 

yj  l^c 

(:>-:^/3^  +  l±^^/3j.(l- |W2|)^^^  (107) 

Finally,  Eq.  (82)  is  no  longer  valid  and  is  replaced  by 

21  V2 

{k  +  2^zf  +  ^  ,  (108) 

for  the  two  remaining  branches.  Growth-rate  curves  for  |u;2l  <  1  will  show  two  unstable 
branches  of  the  dispersion  relation  in  Eq.  (100)  for  large  values  of  real  k. 

If  u»i  7^  0  and  W2  =  0,  then  the  transverse  electromagnetic  components  of  an  eigen- 
mode  are  coupled  through  the  longitudinal  electrostatic  component.  The  large  |^|  behavior 
of  the  dispersion  relation  in  Eq.  (100)  is  similar  to  that  of  the  dispersion  relation  in  Eq.  (69) 
for  the  analogous  case  of  Si  7^  0  and  $2  =  0.  Four  of  the  branches  obey  Eqs.  (79)  and 
(81).  Four  additional  branches  obey  Eqs.  (82)  and  (83)  with  the  quantity  in  Eq.  (84) 
replaced  with  lui,  that  is 


y/2i^c 

(109) 

Cj  ~  k^z  +  1  ± 

■V^Wc 

(110) 
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where 


<=*  =  {t (f‘ -  (1  -  +  [(4 -  (1  -  (I  - 1”'!')]  ^  •  (111) 

The  remaining  two  branches  obey  Eq.  (106).  At  large  real  k,  growth-rate  curves  will  show 
two  unstable  branches  of  the  dispersion  relation  in  Eq,  (100)  if  |u;i|  <  1/V2. 

Finally,  if  neither  wi  nor  W2  vanish,  then  the  large  |fc|behaviors  of  four  of  the  branches 
of  the  dispersion  relation  in  Eq.  (100)  are  given  by  Eqs.  (79)  and  (106).  The  behaviors  of 
the  remaining  six  branches  are  obtained  by  solving  a  cubic  equation  in  (ti)  —  kj3z  —  1)^.  We 
do  not  include  an  analysis  of  the  equation  in  this  paper. 

C,  Numerical  Examples 

Before  the  numerical  examples  are  presented,  one  final  point  of  clarification  must 
be  made.  Throughout  this  paper,  we  have  followed  the  usual  terminology  and  referred 
to  Ei-{k,u))  and  Ei^{k,u})  as  components  representing  RHP-  and  LHP-electromagnetic 
waves,  respectively.  However,  this  nomenclature  is  proper  only  if  Refc  >  0.  If  Refc  <  0,  then 
the  roles  played  by  these  components  are  reversed  and  Ei-{k,Lb)  and  Ei+{k,u))  represent 
LHP  and  RHP  electromagnetic  waves,  respectively.  If  (in  addition)  RetD  >  0,  then  both 
waves  are  backward  traveling.  Moreover,  if  Rew  <  0,  then  both  waves  are  forward  traveling. 

As  in  Sec.  Ill,  the  following  numerical  computations  are  limited  to  the  case  of  real 
k.  Then,  from  Eq.  (102),  Eq.  (100)  is  invariant  under  the  transformation  tD  — >  —u*  and 
k  —>■  —k  —  2/^z.  This  transformation  is  equivalent  to  inverting  a  plot  of  Reu)  vs.  k 
(real)  through  the  point  (fc,Rec2))  =  (-l//3^,0)  and  reflecting  a  plot  of  Im(a;)  vs.  k  (real) 
through  the  vertical  line  k  =  —  1//?^.  Because  of  Eqs.  (103)  and  (104),  no  new  information 
is  obtained  firom  the  transformed  eigenvectors.  Consequently,  the  region  k  <  —lIPz  is 
omitted  from  the  following  plots. 

Parameter  values  used  below  are  the  same  as  those  used  in  previous  numerical  ex¬ 
amples  (i.e.,  Wp/Wg  s=  0.05,  70  =  2,  and  ao  —  0.4).  Functional  forms  considered  for  ’9^(C) 
are  chosen  as  'i'(^)  =  ^(^),  where  ^(^)  is  defined  in  Figs.  3(a)-3(d).'  Moreover  each  of  the 
examples  below  is  the  zmalogue  (for  the  axial-dependent  distribution)  of  the  example  of 
the  same  number  in  Sec.  III.C  (for  the  time-dependent  distribution). 
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Example  1:  For  wi  =  W2  —0,  which  corresponds  to  ’4'(^)  =  $(^)  shown  in  Fig.  3(a), 
the  dispersion  relation  in  Eq.  (100)  reduces  to  the  three  independent  dispersion  relations 

M _ =  0,  M\  I  (k,u)^  =  0,  and  M^z  (k,u>^  =  0  for  uncoupled  right-  and  left- 

circularly  polarized  transverse  waves  and  the  longitudinal  electrostatic  wave,  respectively. 
For  a  given  frequency,  the  respective  wavenumbers  for  these  waves  are  k,  k  +  2fl3z,  and 
ic  +  Ijfiz-  Growth-rate  curves  are  shown  in  Fig.  8(a)  for  =  -1.254  <k  <10.  Plots 
of  Ret2>  vs.  k  for  —l//3z  <  k  <  0  are  presented  in  Fig.  8(b).  Letters  show  corresponding 
points  in  Figs.  8(a)  and  8(b). 


In  Fig.  8(a),  the  growth-rate  curve  segments  OH  and  HI  axe  obtained  from  roots  of 
M++  (k,  ,  so  that  the  corresponding  eigenmodes  have  a  single  nonvanishing  component 

Ei+  (k  +  2/^z,iI^^‘  Referring  to  Fig.  8(b),  we  see  that  Re(2)  >  0  and  k  -|-  2/Pz  >  0  on 
both  segments.  Consequently,  both  segments  represent  unstable  electromagnetic  waves 
that  are  LHP  and  forw2ird  traveling.  All  other  growth-rate  curve  segments  in  Fig.  8(a) 

are  obtained  from  roots  of  M _ ^k,  =  0,  so  that  the  corresponding  eigenmodes  have  a 

single  non  vanishing  component  Bi-  ^k,  u>^ .  For  all  points  of  the  short  growth-rate  curve 
segment  AB,  reference  to  Fig.  8(b)  shows  that  fc  <  0  and  Rea>  <  0.  Consequently,  growth- 
rate  curve  segment  AB  pertedns  to  unstable,  forward-traveling  LHP  electromagnetic  waves. 
Similarily,  k  <  0  and  Red;  >  0  for  eigenmodes  on  growth-rate  curve  segments  BC  and 
DB,  so  that  these  segments  represent  unstable  backward-traveling,  LHP  electromagnetic 
waves.  The  remaining  growth-rate  curve  segments  (BB  and  JK)  pertain  to  unstable 
forward-traveling,  RHP  electromagnetic  waves.  As  expected,  there  is  no  growth  of  the 
RHP  electromagnetic  wave  at  the  resonance  wavenumber  kr  =  1/(1  —  j3z)  =  4.94.  The 
eigenmodes  obtained  from  roots  of  ^k,u)^  are  of  course  completely  decoupled  from 

the  eigenmodes  obtained  from  roots  of  M _ =  0.  Nevertheless,  Fig.  8(a)  represents 

the  limit  approached  by  any  system  with  our  parameters  as  both  Wi  and  W2  approach  zero. 

This  example  is  analogous  to  Example  1  shown  in  Figs.  4(a)  and  4(b)  in  Sec.  III.C. 
Comparing  Figs.  4(a)  and  8(a),  we  see  that  they  differ  in  two  respects.  First,  the  growth- 
rate  curve  in  Fig.  8(a)  obtained  from  M'^  ,  =  0  heis  the  same  form  as  the  growth- 

rate  curve  in  Fig.  4(a)  obtained  from  M++(fc,d))  =  0  but  is  displaced  to  the  left  by 
2/Pz  —  2.508.  Second,  no  information  is  lost  in  Figs.  4(a)  and  4(b)  by  ignoring  the 
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negative  jfe-axis.  However,  only  k  <  -l/fiz  can  be  ignored  in  Figs.  8(a)  and  8(b)  without 
losing  information. 

Example  2:  For  loi  =  0  and  W2  =  2\/2/7r,  which  can  be  obtedned  for  ^(^  =  ^{0 
in  Fig.  3(b)  with  a  =  7r/4,  the  equilibrium  beam  consists  of  two  streams,  eeurh  with  a 
water-bag  distribution  in  ^  of  width  7r/4.  One  distribution  is  centered  at  0  =  mQcz/PzO 
and  the  other  at  ^  =  mQ.ez/pzQ  +  tt.  As  z  varies,  each  center  rotates  about  the  direction 
of  the  applied  field  lines  with  a  characteristic  wavelength  of  ‘I'Kpzo/mQc  =  IttVzq/lOc- 

Growth-rate  curves  (for  —1.25  <  <  16)  and  corresponding  plots  of  Rea)  vs.  k  (for 

-1.25  <  k  <  G)  are  presented  in  Figs.  9(a)  and  9(b).  Plots  of  the  Poynting  ratio  in 
Eq.  (105)  vs  k  (for  -1.25  <  fc  <  16)  appear  in  Fig.  9(c).  Letters  show  corresponding 
points  on  these  plots.  The  letters  have  also  been  chosen  to  correlate  with  letters  on  the 
corresponding  plots  for  Example  2  of  Sec.  III.C  in  Figs.  5(a)-5(c),  which  is  analogous 
to  the  present  example.  Superficially  the  plots  in  Figs.  9(a)-9(c)  are  very  similar  to  the 
corresponding  plots  in  Figs.  5(a)-5(c).  However,  it  is  empheisized  that  the  eigenmodes  are 
very  different  in  the  two  cases.  The  eigenmodes  for  Fig.  5  consist  of  the  nonvanishing 
components  Ei-{k^Cj)  and  Ex+{k,u)  —  2),  whereais  the  eigenmodes  for  Fig.  9  consist  of  the 
nonvanishing  components  Ei-{k,Cj)  2ind  jEi+(fc  +  2/ 

By  comparing  Figs.  9(a)  and  9(b),  it  is  easily  seen  that  the  eigenmodes  belonging  to 
the  growth-rate  curve  segments  BA,  CD,  and  FE  consist  of  LHP,  forward-traveling  elec¬ 
tromagnetic  waves  [from  Ei+{k  +  2/ /3z,u))]  and  LHP,  backward- traveling  electromagnetic 
waves  [from  £?i_(fc,u>)].  [The  corresponding  modes  in  Fig.  5  consist  of  RHP,  forward¬ 
traveling  and  LHP,  backward-traveling  electromagnetic  waves.]  From  Fig.  9(c),  we  see 
that  the  backward  Poynting  flux  is  relatively  strong  for  most  eigenmodes  on  C A  and  that 
it  varies  rapidly  with  k  for  eigenmodes  on  CD  and  FE.  The  infinity  in  the  Poynting  flux 
ratio  at  the  cutoff  at  fc  =  0  is  due  to  the  factor  \{k  +  2Pz)/k\  in  Eq.  (105).  All  of  the 
remaining  growth-rate  curve  segments  in  Fig.  9(a)  pertain  to  eigenmodes  consisting  of  a 
forward-traveling  RHP  cind  a  forward-traveling  LHP  component.  Notice  that  the  branch 
CDH  of  the  dispersion  relation  shows  a  growth  rate  at  the  resonance  kr  —  4.94  which  is 
significantly  greater  than  any  growth  rate  for  the  uncoupled  system  in  Fig.  8(a).  Reference 
to  Fig.  9(c)  shows  that  the  Poynting  flux  of  RHP  electromagnetic  radiation  is  dominant  in 
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the  corresponding  eigenmode.  [This  behavior  is  similar  to  that  foiind  at  k  =  kr  for  Exam¬ 
ple  2  in  Sec.  III.C.]  Figs.  5(a)  and  9(a)  are  very  similar  at  large  values  of  k  in  conformity 
with  previous  analytic  results  pertaining  to  the  large  k  behaviors  Eq.  (69)  when  Si  =  0 
ajid  Eq.  (100)  when  Wi  =  0. 

Example  3:  The  analogue  of  Example  3  of  Sec.  Ill  [whose  stability  properties  are 
summarized  in  Figs.  6(a)-6(d)]  is  obtedned  by  setting  Wi  =  —lij'K  and  W2  —  0  in  Eq.  (98). 
Growth-rate  curves  for  —1.25  <  fc  <  10  are  presented  in  Fig.  10(a).  Deteiils  of  the  growth- 
rate  curves  in  the  negative  k  interval  (—1.25  <  <  0)  are  shown  in  Fig.  10(b).  Corre¬ 

sponding  plots  of  Reu>  vs,  k  (for  —1.25  <  k  <0)  are  presented  in  Fig.  11(a).  Plots  of  the 
Poynting  ratio  in  Eq.  (105)  vs.  k  appear  in  Fig.  11(b).  The  component  Eiz{k  lJPz,u)) 
will  not  necessarily  vanish  for  unstable  eigenmodes  of  this  system.  Consequently,  plots  of 
E\—{k,u})/E\z{k  ■+•  lfPz,Q))  vs.  k  for  unstable  modes  are  presented  in  Fig.  11(c).  Letters 
show  corresponding  points  in  Fig.  10  and  Fig.  11. 

A  detailed  comparison  of  Figs.  10  and  11(a)  gives  the  following  description  of  the  un¬ 
stable  eigenmodes.  The  components  of  an  eigenmode  pertaining  to  the  short  growth-rate 
curve  segment  DR  are  two  forward-traveling  LHP  electromagnetic  waves  [from\Ei_(fc,a)) 
and  £?i+(fc-f  w)]  and  a  backward-traveling  electrostatic  wave  [from  Eu[k+\J ^z^Cj)]. 

The  components  pertaining  to  the  growth-rate  curve  segments  REF,  ABCG,  and  MN 
are  a  backward-traveling  LHP  electromagnetic  wave  [from  Ei-{k,u))],  a  forward-traveling 
LHP  electromagnetic  wave  [from  Ei+{k  -f  2//3z,£j)],  and  a  forward  traveling  electrostatic 
wave  [from  Eiz{k  +  l//?^,a))].  All  other  growth-rate  curve  segments  have  eigenmodes  con¬ 
sisting  of  forward  traveling  LHP  and  RHP  electromagnetic  waves  and  a  forward  traveling 
electrostatic  wave. 

The  most  rapidly  growing  eigenmode  of  this  system  is  that  at  point  C  in  Figs.  10(a) 
and  10(b).  The  components  of  this  eigenmode  axe  a  backward-traveling  LHP  wave,  a 
forward-traveling  LHP  wave,  and  a  forward-traveling  electrostatic  wave.  Reference  to 
Figs.  11(b)  and  11(c)  shows  that  the  backward-traveling  component  [which  arises  from 
Ei-{k,u))]  is  the  largest  component  both  in  amplitude  and  energy  treinsfer.  A  similar 
situation  was  found  for  small  |fcl  in  Fig.  6(a)  for  the  axial-dependent  case  except  that  the 
forward-traveling  electromagnetic  component  was  found  to  be  RHP. 
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Two  branches  of  the  growth-rate  curves  in  Figs.  10(a)  and  10(b)  {MNO  and 
ABCGHI)  show  moderate  growth  rates  at  the  resonance  wavenumber  kr  =  4.94.  For 
both  of  these  branches,  Figs.  11(b)  and  11(c)  show  that  the  eigenmode  at  ^  =  4.94  has  a 
relatively  large  RHP  electromagnetic  component.  Again  we  see  that  this  behavior  differs 
from  that  of  the  xincoupled  system  in  Fig.  8(a) ,  which  shows  no  growth  of  RHP  electromag¬ 
netic  waves  at  the  resonance  wavenumber.  The  growth  rates  at  kr  are  approximately  the 
same  in  Fig.  6(a)  for  the  time-dependent  equilibrium  and  Fig.  10(a)  for  the  axial-dependent 
equilibrium;  however,  the  electrostatic  components  of  the  corresponding  eigenvectors  are 
of  greater  relative  amplitude  in  the  axial-dependent  case  than  in  the  time-dependent  case. 
[Compare  Fig.  11(c)  with  Fig.  6(d).] 

Finally,  at  large  values  of  k,  Figs.  10(a)  and  6(a)  approximate  each  other  closely.  This 
fact  conforms  with  our  previous  results  giving  the  large-|fc|  behaviors  of  Eq.  (69)  for  S2  =  0 
and  Eq.  (100)  for  W2=0-  However,  the  corresponding  eigenmodes  [E  in  Eq.  (60)  and  E'  in 
Eq.  (91)]  are  different  even  in  the  limit  of  large  k.  By  comparing  Fig.  6(d)  with  Fig.  11(c), 
it  is  seen  that  (at  large  k)  the  electrostatic  component  is  relatively  much  stronger  in  the 
ceise  of  the  axial-dependent  equilibrium  distribution. 

Example  4'  To  obtaiu  the  analogue  of  Example  4  of  Sec.  III.C  (whose  stability  prop¬ 
erties  are  summarized  in  Fig.  7),  choose  wi  —  —  1.  These  values  are  obtained  by 

choosing  ^'(C)  =  IZ-oo^(C  “  2n7r),  so  that  in  effect  the  phase  of  any  particle  is  given 
by  (^  =  rnJ^lczJpz).  Growth-rate  curves  (for  -1.254  <  k  <  8)  and  corresponding  plots 
of  Red;  vs.  k  (for  —1.254  <  ^  <  0)  are  presented  in  Figs.  12(a)  and  12(b),  respectively. 
For  unstable  branches  of  the  dispersion  relation  in  Eq.  (100),  plots  of  the  Poynting  flux 
ratio  in  Eq.  (105)  vs.  k  and  £i-(fc,w)/£^iz(fc  +  l//3*,w)  vs.  k  are  presented  in  Figs.  12(c) 
amd  12(d),  respectively.  Letters  show  corresponding  points  in  these  plots.  The  letters 
correspond  only  loosely  to  those  in  Fig.  7. 

Comparing  Figs.  12(a)  and  12(b),  we  see  that  the  components  of  an  eigenmode 
on  growth-rate  curve  segment  DN  are  a  backward-traveling  LHP  electromagnetic  wave 
[Ei+(^  +  2/(3z,<^)],  a  forward-traveling  LHP  electromagnetic  wave  [Ei_(fc,d))],  and  a 
backward-traveling  electrostatic  wave  [Eiz{k  l//?^,^))].  Eigenmodes  on  growth-rate 
curve  segments  ACB  and  NEF  consist  of  a  backward-traveling  LHP  electromagnetic 
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wave  [^i_(fc,a>)],  a  forward- traveling  LHP  electromagnetic  wave  [Ex+{k-^2/ and  a 
forweird-traveling  electrostatic  wave  \Exz{k  l//3a,a>)].  Eigenmodes  on  all  other  growth- 
rate  curve  segments  consist  of  forward  traveling  RHP  and  LHP  electromagnetic  waves  and 
a  forward-traveling  electrostatic  wave. 

Some  properties  of  the  growth-rate  curves  in  Fig.  12(a)  are  similar  to  those  in  Fig.  7(a). 
T-ilfft  Fig.  7(a),  Fig.  12(a)  shows  no  growth  at  large  values  of  k.  Both  sets  of  curves  show 
very  large  growth  rates  at  small  values  of  |fc|,  where  backward  waves  occur  [i.e.,  near  point 
B  in  Fig.  7(a)  and  point  C  in  Fig.  12(a)].  The  eigenmode  at  point  C  in  Fig.  12(a)  consists 
of  a  backward- traveling  LHP  electromagnetic  wave  [from  Ei_(fc,(D)],  and  forward-traveling 
LHP  electromagnetic  and  electrostatic  modes.  Reference  to  Figs.  12(c)  and  12(d)  shows 
that  the  backward-traveling  component  exceeds  the  other  two  components  in  amplitude. 
In  Sec.  III.C,  a  similar  situation  was  found  to  exist  at  point  B  in  Fig.  7(a),  except  that 
the  forward-traveling  electromagnetic  component  is  RHP.  Like  Fig.  7(a),  Fig.  12(a)  shows 
a  fairly  large  growth  rate  at  the  resonance  wavenumber  on  the  branch  ACBIJ.  Moreover, 
Figs.  12(c)  and  12(d)  show  that  the  corresponding  eigenmode  has  a  relatively  strong  RHS 
electromagnetic  component. 

Finally,  notice  that  no  gap  appears  in  the  growth-rate  curves  in  Fig.  12(a)  to  corre¬ 
spond  to  the  gap  between  points  F  and  H  in  Fig.  7(a). 

D.  Remarks  Concerning  Numerical  Examples. 

The  analysis  of  the  above  numerical  examples  for  the  axial-dependent  equilibrium 
leads  to  the  same  general  conclusions  as  those  given  in  Sec.  III.C  for  the  time-dependent 
equilibrium.  A  rich  structure  of  different  growth-rate  curves  and  unstable  eigenmodes 
can  be  induced  by  varying  the  form  of  'i'(C))  the  values  of  wi  and  W2.  A  suitable 
choice  of  '5'(C)  can  significantly  increase  growth  rates  of  RHP  electromagnetic  waves  at 
the  resonance  wavenumber  k  —  1/(1  —  Pz)  aJid  can  significantly  increase  or  reduce  growth 
rates  at  large  values  of  k.  At  small  values  of  k,  where  backward-traveling  components  are 
present,  growth  rates  and  the  structures  of  eigenmodes  depend  strongly  on  the  form  of 
’4^(C).  This  fact  suggests  that  properties  of  absolute  instabilities  may  depend  strongly  on 
5'(C).  However,  a  study  of  this  conjecture  has  not  been  carried  out. 
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For  the  same  parameters  (^^c,7o,<»o),  growth  rate  curves  for  corresponding  (si  =  Wu 
S2  =  1V2)  time-dependent  and  axial-dependent  systems  usually  show  some  resemblance. 
Nevertheless,  the  eigenmode  structures  are  very  different  in  the  two  cases.  In  the  time- 
dependent  case,  unstable  eigenmodes  for  coupled  systems  consist  of  two  or  three  compo¬ 
nents  of  the  same  wavenumbers  and  different  frequencies,  whereas  in  the  axial-dependent 
case  the  components  have  the  same  frequencies  and  different  wavenumbers.  Moreover, 
for  small  values  of  k,  the  handedness  and  directions  of  motion  of  components  may  differ 
between  the  two  cases.  . 
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V.  Conclusions 


We  have  studied  stability  properties  of  a  relativistic  electron  beam  propagating  along 
an  applied  magnetic  field  using  the  Maxwell- Vlasov  equations  under  the  constraint 

that  spatially  dependent  quantities  are  functions  of  z  only.  Of  particular  interest  are 
those  cases  in  which  the  equilibrium  distribution  is  not  uniformly  random  in  the  electron 
gyro-phase  angle  Two  equilibrium  distributions  have  been  considered.  These  are  the 
time-dependent  distribution  fo(p±iPz,0>  where  ^  =  <f>  —  Clct/lf}  and  the  spatial-dependent 
distribution  /o(p±)P«>C)>  where  C  =  ^  “  niQczIPz-  Since  neither  of  these  distributions 
can  be  converted  into  the  other  by  a  Lorentz  transformation,  the  distributions  represent 
two  physically  different  systems.  It  is  found  that  in  general  the  Fourier  components  of  the 
perturbed  electric  and  magnetic  fields  are  related  the  integral  equations  (26)-(28)  for  the 
case  of  the  time-dependent  equilibrium  distribution,  and  by  the  integral  equations  (41)-(43) 
for  the  case  of  the  spatial-dependent  eqmlibrium  distribution.  In  our  numerical  analysis, 
however,  we  consider  special  cases  in  which  the  integral  equations  reduce  to  algebraic 
equations  even  though  the  equilibrium  distribution  is  not  uniformly  random  in  phase. 

If  there  is  no  spread  in  electron  energies  (or  equivalently  p)  in  the  time-dependent 
equilibrium  distribution,  then  the  integral  equations  (26)-(28)  reduce  to  just  three  algebraic 
equations  [Eq.  (58)]  relating  the  Fourier  components  Ei-{k,Lo),  —  2uJc),  and 

Eiz{k,u>  —  ojc)  of  tlio  perturbed  fields.  Consequently,  an  eigenmode  of  the  system  consists 
of  a  RHP  electromagnetic  wave,  a  LHP  electromagnetic  wave,  and  an  electrostatic  wave. 
These  components  have  the  same  wavenumber,  and  the  same  spatial  and  temporal  growth 
or  decay  rates,  but  have  different  frequencies.  [The  electrostatic  component  is  decoupled 
if  the  Fourier  coefficient  giipo,<^)  in  Eq.  (57)  vanishes,  and  all  three  components  decouple 
if  52  (Po>  o)  also  vanishes.] 

If  there  is  no  spread  in  the  axial  component  of  momentum  (pz)  in  the  spatial-dependent 
equilibrium  distribution,  then  the  integral  equations  (41)-(43)  reduce  to  just  three  al¬ 
gebraic  equations  [Eq.  (89)]  relating  the  perturbed  field  Fourier  components  Ei-{k,u), 
Ei+{k  +  2mClc/PzO,<*’)t  Eiz{k  +  mQ,c/PzO,<^)-  Therefore,  the  components  of  an  eigenmode 
are  a  RHP  electromagnetic  wave,  a  LHP  electromagnetic  wave,  and  an  electrostatic  wave. 
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These  components  have  the  same  frequency,  and  the  same  spatial  and  temporal  growth  or 
decay  rates,  but  have  different  wavenumbers.  [In  analogy  with  the  time-dependent  case, 
the  electrostatic  component  is  decoupled  if  the  Fourier  coefficient  hi(pj.,Pzo)  in  Eq.  (88) 
vanishes,  and  all  of  the  components  decouple  if  ^(p±,Pzo)  also  vanishes.] 

Numerical  computations  of  stability  properties  have  been  carried  out  for  both  the  time- 
and  spatial-dependent  equilibrium  distributions  for  the  case  where  no  spread  is  present 
in  both  p  and  the  pitch  angle  o  (or  equivalently  in  both  p±,  and  pz)-  In  this  case  the 
frequencies  and  wavenumbers  can  be  normalized  to  the  relativistic  cyclotron  frequency 
by  HpfiniTig  d>  =  (jj/uc  and  k  =  ckjuj^.  The  computations  are  restricted  to  real  values  of  fc, 
so  that  InuD  >  0  indicates  temporal  growth.  It  is  found  that  (for  fixed  applied  magnetic 
field,  energy,  and  pitch  angle)  a  rich  variety  of  growth-rate  curves  and  eigenmodes  can  be 
obtained  by  changing  the  dependence  of  the  equilibrium  distribution  on  the  phase  angle. 
Appropriate  choices  of  the  phase-angle  dependence  can  significantly  increase  growth  rates 
near  the  resonance  wave  number  kf.  =  1/(1  “  i^^).  Growth  rates  at  large  values  of  k  can 
on  the  one  hand  be  enhanced  and  on  the  other  hand  be  suppressed  altogether.  Moreover, 
finite  intervals  (in  k)  of  no  growth  can  be  produced.  Finally,  growth  rate  curves  and 
the  form  of  eigenvectors  at  small  values  of  |fc|,  where  backward  traveling  components  are 
present,  are  particularly  sensitive  to  the  (/^dependence  of  the  equilibrium  distribution. 

As  an  important  area  in  our  current  research,  we  are  analyzing  the  integral  equations 
to  determine  the  structures  of  eigenmodes  in  the  general  case. 
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C  PROGRAM  TO  COMPUTE  CORRELATION  DIMENSION  FROM  DATA  SET 
C 

C  INPUT  MATRIX  (ORBIT  TRAJECTORY  )  FROM  FILE 
C  INPUT  IDELA Y  THE  TIME  DELAY  IN  THE  CONSTRUCTION  OF 
PSUEDO-VECTOR 

C  INPUT  IMBED  THE  DIMENSION  OF  THE  PSUEDO-VECTOR 

C  INPUT  NUM  THE  NUMBER  OF  POINTS  FOR  WHICH  THE  CORRELATION 

INTEGRAL 

C  WILL  BE  PLOTED 

REAL  A,C,RJ)MAX,DMIN 
DIMENSION  A(5000),C(500),R(500) 

INTEGER  IMBED,IDELAY,NUM 
NAMELIST/INPUT/IMBED,IDELAY,NUM 
OPEN  (2,FILE=’C-INPUT.DAr,STATUS=’OLD’) 

READ  (2,INPUT,ERR=1000) 

WRITE  (6,INPUT) 

OPEN  (4T^ILE=PLOT.DAr,STATUS='OLD’) 

OPEN  (9  JILE='CPLOT.DAr,STATUS='NEW’) 

DO  12  1=1,5000 

READ  (4,10, IOSTAT=IO)  A(I) 

IF  (10  .EQ.  -1)THEN 
NDATA  =  I-1 
GO  TO  13 
ENDF 

12  CONTINUE 

10  FORMAT  (:,(E15.6)) 

C  FIND  THE  NUMBER  OF  DATA  POINTS  IN  THE  INPUT  FILE 
C  NDATA  =  0 
C  DO  11  1=1,5000 
C  IF  ( A(I)  .NE.  0.0)THEN 
C  NDATA  =  NDATA+ 1 
C  ELSE 
C  GO  TO  13 
C  ENDIF 

11  CONTINUE 

13  CONTINUE 
C 

C  FIND  THE  MAXIMUM  AND  MINIMUM  DISTANCE  BETWEEN  ELEMENTS  IN 
THE  DATA-SET 
C 

CALL  MAXMIN(A,NDATA,DMAX,DMIN,IMBED,IDELAY) 
WRITE(6,*)'NDATAMINANDMAXVALUES',NDATA,'  ’,DMIN,'  ',DMAX 

C  NOW  FIND  THE  CORRELATION  INTEGRAL 
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CALLCORREL(A,NDATAJDMAX,DMIN,IMBED,IDELAY,NUM,C,R) 

C  NOW  WRITE  THE  CORRELATION  DATA  TO  AN  OUTPUT  FILE 
DO  22  I=1,NUM 

WRITE(9,20)R(I),C(I)ALOG(R(I)),ALOG(C(I)) 

WRITE(6,20)  R(I),C(I),ALOG(R(I)),ALOG(C(I)) 

20  FORMAT(4(E15.6)) 

22  CONTINUE 

STOP  TsfORMAL  TERMINATION* 

1000  STOP  *ERROR  IN  NAMELIST  READ' 

END 

SUBROUTINE  MAXMIN(A^ENGTH,DMAXT)MIN,IMBED,IDELAY) 

REAL  A(NLENGTH),DMAX,DMIN,Z,D,D1 
INTEGER  IMBED,IDELAY 

C  THIS  IS  A  STATEMENT  FUNCTION 
Z(I,J)  =  A(I+(J-1  )*IDELAY) 

NC  =  NLENGTH  -  (IMBED-1)*IDELAY 
DMAX  =  0.0 
DMIN=I.0E3 
DO30I=l,NC-l 
DO20K=I+l,NC 
D  =  0.0 

DO  10  JJ=1,IMBED 
D1  =(Z(I,JJ)-Z(K,JJ))**2 
D  =  D1+D 
10  CONTINUE 

DMAX  =  MAX1(DMAX,D) 

DMIN  =  AMIN1(DMIN,D) 

20  CONTINUE 
30  CONTINUE 

DMAX  =  SQRT(DMAX) 

DMIN  =  SQRT(DMIN) 

RETURN 

END 

SUBROUTINE  CORREL(A,NDATA,DMAX,DMIN,IMBED,IDELAY,NUM,C,R) 
REAL  A(NDATA)TDMAX,DMIN,DINC,C(NUM),R(NUM) 

INTEGER  IMBED,IDELAY,NUM 

C  THIS  IS  A  STATEMENT  FUNCTION 
Z(I,J)  =  A(I+(J-1)*IDELAY) 

DINC  =  (DMAX  -  DMIN)/FLOAT(NUM-1.0) 

DO  10  I=1,NUM 
C(I)  =  0.0 
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R(I)  =  DMIN  +  (I-1)*DINC 
10  CONTINUE 

NC  =  NDATA  -  (IMBED-1)*IDELAY 
DO  40  I=1,NC-1 
DO30K=I+l,NC 
0=0.0 

DO  20  JJ=1, IMBED 
D1=(Z(I,JJ)-Z(K,JJ))**2 
D  =  D  +  D1 
20  CONTINUE 

MM  =  INT(  (SQRT(D)  -  DMIN)/  DINC ) 
DO  25  M1=MM+1,NUM 
C(M1)=1+C(M1) 

25  CONTINUE 

30  CONTINUE 
40  CONTINUE 
DO50L=l,NUM 
C(L)  =  C(L)/NC**2 
50  CONTINUE 
RETURN 
END 
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C  Program  ssvdly.for 

c  PROGRAM  TO  EVALUATE  LYPUNO V  SPECTRUM 

C  FROM  TRAJECTORY  DATA 

c 

c  The  trajectory  data  in  this  revision,  comes  from 
c  the  psuedo-trajectories  generated  by  the  svd 
c  of  a  single  observable  time  series.  svdplot.for 
c 

c  MODIFICATION  OF  CODE  BY  WOLFE  .ET  AL 

c  These  final  modification  are  implemented  on  this 
c  the  16  day  of  March  1995,  and  substainially  enhances 
c  and  modifies  the  results  of  Wolfe  .et.  al 

C 

INTEGER  DIM,EVOLV,IND,BACKIT,BACKMAX,NPT 

REAL  PT,VOUT,VPREV,VOL 

REAL  VOLNEW,NORM,duml,dum2 

REAL  Z,SUM,ZLYP,test 

REALDOT,dum3T)EX 

Dimension  PT(5,5),Vout(5),yprev(5,5),Vol(5),IND(6) 
Dimension  Volnew(5)T)uml(5000),dum2(5) 

Dimension  Z(5000,5),Sum(5),Zlyp(5),DEX(5000),Dum3(5) 
LOGICAL  GOODPT 
COMMON/BLKl/  Z,IND,DEX,DIM 
C  TfflS  PERMITS  EMBEDING  DIMENSION  UP  TO  5 
NAMELIST/NAMEIN/  DT,SCALMN,EVOLV,BACKMAX 


zeroeff=  l.Oe-5 
C 

C 

OPEN  ( 1  ,FILE=’svdata.dat',STATUS='OLD') 

OPEN(2,FILE='sdlyp_in.dat',STATUS='OLD') 

OPEN(5,FILE=’lypdspec.dat',STATUS='NEW) 


READ  (2,NAMEIN,ERR=4000) 

WRITE(6,NAMEIN) 

C  ' 

C  IND(l)  POINTS  TO  THE  FIDUCIARY  TRAJECTORY 
C  IND(2)  POINTS  TO  THE  SECOND  TRAJECTORY,  etc 
C  SUM  HOLDS  THE  RUNNING  EXPONENT  ESTIMATE  SANS  1/  TIME 

C  ITS  IS  THE  TOTAL  NUMBER  OF  PROPAGATION  STEPS 

IND(1)=1 


AII-4 


o  n  n 


APPENDIX  n 


BACKIT  =  0 
DO  2  J=1,DIM 
SUM(J)  =  0.0 
2  CONTINUE 

ITS  =  0 

READ  IN  THE  TIME  SERIES 

Read(l,*)NPT,DIM 
D04I=1,NPT 

READ(1,*)(Z(I,K),K=1,DIM) 

4  CONTINUE 
C 

C  FIND  THE  NEAREST  NEIGHBOR  TO  THE  FIRST  DATA  POINT 

5  D0  1000JDIM=1,DIM 

DO20I=l,NPT 
DO  =  0.0 

DO10J=l,DIM 

DO  =  DO +(Z(I,J)  -  Z(IND(JDIM),J))**2 
10  CONTINUE 

duml(i)  =  D0 
DEX(i)  =  i 

Write(*,*)  DEX(i)  and  i',Dex(i),i 
20  CONTINUE 
c  do  21  iii=l,npt 
c  duml(iii)  =  d(jdim,iii) 

c  DEX(ni)  =  FLOAT(m) 

c21  continue 

c  CALL  SORT(dum  1  ,INDEX,SCALMN,NPT) 

CALL  SORT2(NPT,DUMl,DEX) 

C  THIS  WILL  RESULT  IN  AN  INDEX  OF  POINTS  OUTSIDE  OF  THE  NOISE 

C  SCALE  LENGTH  AND  IN  ORDER  OF  INCREASING  DISTANCE  FROM  THE 

C  FIDUCIARY  POINT 

C  DETERMINE  IF  A  SELECTED  POINT  IS  A  GOOD  POINT,  I.E.  IT  IS 
C  SEPARATED  BY  AT  LEAST  THE  SCALE  MINIMUM  FROM  THE  ALL 

PREVIOUSLY 
C  SELECTED  POINTS 

INEXT=1 

30  DO  50  I=INEXT,NPT 
GOODPT  =  .TRUE. 
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DO40J=l,JDIM 
Call  DIST(J,IJ)0) 

GOODPT  =  DO  .GT.  SCALMN  .AND.  GOODPT 
40  CONTINUE 
IF(GOODPT)THEN 
INDEX  =  IFIX(DEX(I)) 

IND(JDIM+1)  =  INDEX 
INEXT  =  I+1 
GO  TO  60 
ENDff 

50  CONTINUE 
C 

C  TfflS  LINE  OF  CODE  IS  REACHED  ONLY  IF  A  CHECK  OF  ALL  NPTS 

YIELD 

C  NO  GOOD  POINTS 

STOP  NOT  ENOUGH  DATA’ 

60  CONTINUE 
C 

C  NOW  CHECK  THAT  THE  POINT  SELECTED  IS  NOT  DEGENERATED 

C  GIVE  THE  LOCAL  NAME  OF  THE  VECTOR  VOUT 

DO  70  n=l,DIM 

VOUT(II)  =  Z(IND(JDIM+l),n)  -  Z(IND(1),II) 

70  CONTINUE 

CALL  GRAM(  VOUT,  VPREV,DIM,JDIM-1, test) 

NORM  =  test 

IF(NORM  .LT.  zeroeff  )THEN 

F(BACKIT  .GE.  BACKMAX)  STOP  IVIAX  BACK  SUBSTITUTION 
REACHED' 

BACKIT  =  BACKIT+  1 

GO  TO  30 

ENDF 

C  THIS  LINE  IS  REACHED  ONLY  ff  NON-DEGENERATE 

ORTHOGONALIZATION 

C  HAS  TAKEN  PLACE,  AS  IS  DESIRED 

VOL(l)  =  NORM 
ff(JDIM  .EQ.  1)  GO  TO  1000 

C  STORE  THE  COMPONENTS  OF  VPREV  IN  A  DUMMY  ARRAY  (DUM2) 
C  INORDER  TO  PASS  IS  TO  THE  SUBROUTINE 

D0  7ini=l,DIM 
DUM2(m)  =  VPREV(JDIM,ni) 

DUM3(in)  =  VPREV(l,in) 
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71  CONTINUE 

CALL  VDOT(VOUTJ)UM2,DlM, DOT) 
c  VOL(JDIM)  =  VOL(JDIM-l)*ABS(DOT) 

VOL(JDIM)  =  ABS(DOT) 

c  Compare  the  results  of  NORM  and  DOT,  they  should  have  the  same  values 

c  Write(*,*)  'Jdim  and  NORM  and  DOT, JDIM,  NORM(JDIM),  DOT 

c  Check  the  orthogonalization  of  vectors 

c  CALL  VDOT(Dum3,Dum2,Dim,dot) 
c  Write(*,*)  'The  dot  Product  of  e-1  and  e-JDIM  is  ',dot 

C 

1000  CONTINUE 


C  THE  COORDINATES  OF  THE  EVOLVED  POINTS  ARE  AS  FOLLOWS 


DO  1010J=14)IM+1 

F(IND(J)+EVOLV  .GT.  NPT)STOP  'END  OF  FILE  REACHED' 

DO  1020  I=1,DIM 
PT(J,I)  =  Z(IND(J)+EVOLV,I) 

1020  CONTINUE 
1010  CONTINUE 
C 

C  FIND  THE  VOLUMES  OF  THE  EVOLVED  POINTS 

DO  2000  JDIM=1,DIM 
DO  2010  II=1,DIM 
VOUT(n)  =  PT(JDIM+l,n)  -  PT(1,II) 

2010  CONTINUE 

CALL  GRAM(VOUT,VPREV,DIM,JDIM-l,TEST) 

NORM  =  TEST 
VOLNEW(l)  =  NORM 
IF(JDIM  .EQ.  1)  GO  TO  2000 

C  STORE  THE  COMPONENTS  OF  VPREV  IN  A  DUMMY  ARRAY  (DUM2) 

C  INORDER  TO  PASS  IS  TO  THE  SUBROUTINE 

DO  171  m=l,DIM 
DUM2(ni)  =  VPREV(JDIM,in) 

171  CONTINUE 

CALL  VDOT(VOUTT)UM2,DIM,DOT) 
c  VOLNEW(JDIM)  =  VOLNEW(JDIM-l)*ABS(DOT) 
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VOLNEW(JDIM)  =  ABS(DOT) 
c  Compare  the  NORM  and  DOT  for  the  New  Volumes 

c  Write(*,*)  TJew  Vols:  JDIM  and  NORM  and  DOr,JDIM,NORM(JDIM)JDOT 
2000  CONTINUE 

C  EVALUATE  APPROXIMATE  EXPONENTS 
ITS  =  ITS+1 
DO  3000  J=1,DIM 

SUM(J)  =  SUM(J)  +  ALOG(VOLNEW(J)AVOL(J))/(FLOAT(EVOLV) 

.  *DT*ALOG(2.0eO)) 

ZLYP(J)  =  SUM(J)/FLOAT(ITS) 

3000  CONTINUE 

WRITE(5,*)  (ZLYP(J),J=1,DIM) 

55  format(il0,2x,3(el6.5)) 

C 

C  FIND  THE  DIM-NEAREST  'GOOD'  NEIGHBORS  OF  THE  EVOLVED 
C  FIDUCIARY  POINT 

IND(  1 )  =  IND(  1 )  +  EVOLV 
IF(  IND(l)  XT.  NPT)THEN 
GO  TO  5 
ENDIF 

STOP  'NORMAL  EXECUTION 
4000  STOP  ERROR  IN  NAMELIST 
END 

SUBROUTINE  DIST(I,J,D0) 

INTEGER  IND,INDEX  J)IM 
REALZ,D0,DEX 

Dimension  IND(6)^(5000,5),DEX(5000) 

COMMON/BLKl/ Z,IND.DEX,DIM 

DO  =  O.ODO 

INDEX  =  IFD((DEX(J)) 
do  5  ii=l,DIM 

DO  =  DO  +  ( Z(IND(I),ii)  -  Z(INDEX,ii)  )**2 
5  continue 

DO  =  SQRT(D0) 

END 


SUBROUTINE  SORT(A,INDEX,SCALMN,NLIST) 
INTEGER  NLIST,INDEX,I,IMIN 
REALA,DUM,D 
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REAL  SCALMN 

DIMENSION  INDEX(5000),A(5000),Dum(5000) 
c  intialize  the  dummy  variable  dum 
do  10  ii=l,NLIST 
DUM(ii)  =  A(ii) 

10  continue 

DO  100  K=1,NLIST 
D=1.0E10 
DO  50  I=K,NLIST 

IF(DUM(I)  .LE.  D  .AND.  DUM(I)  .GE.  SCALMN)THEN 
D  =  DUM(I) 

IMIN  =  I 
ENDIF 

50  CONTINUE 

DUM(IMIN)  =  DUM(K) 

DUM(K)  =  D 
INDEX(K)  =  IMIN 
100  CONTINUE 
RETURN 
END 

SUBROUTINE  GRAM(VIN,VPREV,DIM,NPREV,NORM) 

INTEGER  DIM,NPREV 

REAL  VIN,DUM,GSC,VPREV,NORM,DUM2,test 
Dimension  Vin(5),Dum(5),Gsc(5),Vprev(5,5)dDum2(5) 

C 

C  INITIALIZE  DUM 

DO  10  I=1,DIM 
DUM(I)  =  VIN(I) 

10  CONTINUE 

C 

C  FIND  THE  COEFFICIENTS 

IF(NPREV  .EQ.  0)  GO  TO  40 
DO  20  J=1,NPREV 

C  STORE  THE  COMPONENTS  OF  VPREV  IN  A  DUMMY  ARRAY  (DUM2) 

C  mORDER  TO  PASS  IS  TO  THE  SUBROUTINE 

DO  71  ni=l,DIM 
DUM2(ni)  =  VPREV(J,in) 

71  CONTINUE 

CALL  VDOT(VIN.DUM2,DIM,test) 

GSC(J)  =  test 
20  CONTINUE 

DO  30  L=1,NPREV 
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D0  30J=1J)IM 

DUM(J)  -  DUM(J)  -  GSC(L)*VPREV(L,J) 
30  CONTINUE 

40  CALL  VD0T(DUMJ)UM,DIM,N0RM) 
IF(NORM  .EQ.  0)  RETURN 
NORM  =  SQRT(NORM) 

DO  50  J=1,DIM 

VPREV(NPREV+1  ,J)  =  DUM(J)/NORM 
50  CONTINUE 
RETURN 
END 

SUBROUTINE  VD0T(V1,V2JDIMJ30T) 
INTEGER  DIM 
REAL  Vl,V2J)OT 
Dimension  V1(5),V2(5) 

DOT  =  0.0 

DO  10  I=1,DIM 

DOT  =  DOT  +  Vl(I)*V2(I) 

10  CONTINUE 
RETURN 
END 
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